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Mass Transfer Through Compressible 


Turbulent see Layers 


T. K. SHERWOOD? and H. 


New data are presented on mass transfer to air at velocities up to 


520 m. per sec. (Mach 2.02). 


These were obtained by measuring the 


rate of sublimation of naphthalene, thymol, and camphor from cylinders 
placed longitudinally in air streams in three subsonic and supersonic 


ducts. 


Earlier published data are limited to air velocities of about 10 m. 


per sec. 
following Chilton and Colburn. 
this correlation. 


These are collected and correlated as a graph of jp vs. Re, 
The new subsonic data are found to fit 


At supersonic air velocities there is an important effect of Mach 


number, superimposed on the effect of Reynolds number. 


Both sub- 


sonic and supersonic data are correlated by introducing a function of 
Mach number similar to that found to apply to skin friction at supersonic 


velocities. 


The results lend semi-quantitative support to the conclusion 


that the effect of Mach number is the same for mass transfer and for 


skin friction. 


HE nature of turbulence and of 

flow past surfaces is now suffici- 
ently well understood that useful rela- 
tions between mass transfer, heat 
transfer, and skin friction may be 
obtained from the basic theory. These 
three phenomena, once considered 
separate fields for empirical study, 
are now recognized as being related 
by their common dependence on the 
mechanics of flow. For engineering 
use the relations between heat, mass, 
and momentum transfer are expressed 
in the form of modifications of the 
Reynolds analogy, as developed by 
Chilton and Colburn, Von Karman, 
and others. 

Since World War II the increased 
interest in high speed flight has 
stimulated both theoretical and ex- 
perimental studies of the effect of gas 
compressibility on skin friction. 
Studies of heat transfer at supersonic 
velocities have also been carried out, 
and the two processes related. Very 


little has been done on mass transfer 
at high air speeds, though the experi- 
mental study of mass transfer is 
potentially simpler than either heat 
transfer or friction tests as a method 
of investigating compressible flow 
phenomena. 

The present paper describes mea- 
surements of the rate of sublimation 
of a solid from a cone-cylinder test 
piece in air at velocities up to 520 m. 
per sec. (Mach 2.02) and discusses 
the relation of the results to published 
skin friction measurements. 


The Relations between Mass, 
Momentum, and Heat Transfer 
at Subsonic Velocities 

All three processes depend on both 
molecular motion and the mixing of 
eddies being interchanged by turbu- 
lent streams. Near a solid surface the 
eddy motion is weak or absent, and 
transfer depends on molecular motion. 
In the main turbulent stream eddy 
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transfer is rapid and molecular motion 
contributes little. The overall process 
depends both on the molecular pro- 
perties and on the conditions of tur- 
bulence and flow. 

The Reynolds “analogy”, relating 
friction and heat transfer, is based on 
the similarity of the mechanisms of 
heat and momentum transfer by 
eddies, and makes no allowance for 
the role of molecular motion. Later 
modifications of the Reynolds rela- 
tion, some of which are summarized 
in an earlier paper‘!), provide for the 
decreasing importance of eddy trans- 
fer and the dominant role of molecular 
transport near the wall, using the 
generalized Nikuradse-Von Karman 
velocity-distance correlation to ob- 
tain eddy transfer rates. As developed 
for mass transfer, these are all of the 


form 
ke 
=f E : é (1) 
CrVo pD E 


where k, is the mass transfer coeffici- 
ent based on a concentration driving 
force, C, is the local skin friction 
coefficient, and V, the free stream 
velocity. In gases the potential is or- 
dinarily expressed in terms of partial 
pressure of the diffusing substance, 


and k, is replaced by k,RT. 


The first of the dimensionless groups 
in the function on the right is the 
Schmidt group, which will be repre- 
sented by Sc. The second is the ratio 
of the eddy “viscosity”, «, to the mass 
eddy diffusivity; this ratio appears to 
be a constant equal to about 0.70. 
The various theories lead to different 
forms of the function on the right. 
The simplest of these, and one which 
has found wide acceptance, is Sc?/8, 
due to Chilton and Colburn. 

The corresponding relation for heat 
transfer involves the heat transfer 
coefficient and the Prandtl group. 
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Combining the two, there results 
kcRT 


Sc?2/3 ( jp) = 
Vo P ‘ 
h 2 C; 
— Pr?/3 (= jy) = (2 
CpVop f 2 
The three dimensional groups of 


Equation 2 have been shown to be 
related to the Reynolds number in 
approximately the same way for a 
variety of cases of turbulent flow of 
gases and liquids past solid surfaces‘). 
The first is commonly represented by 
the symbol j, and the second by jy. 


Evaporation from Flat Plates 

Data on evaporation from flat plates 
have been collected from the liter- 
ature and plotted in Figure 1, not 
only as an illustration of the applic- 
ability of Equation 2, but for com- 
parison with the rates of sublimation 
at high air speeds as measured in the 
work to be described. 

The data shown are plotted as j,, 
vs. Rep, where Re, is the Reynolds 
number based on free-stream condi- 
tions and the down-wind length of 
the wetted surface. Most of the data 
were obtained by placing flat pans of 
liquid flush with the floor of a wind- 
tunnel, though Maisel, Powell, and 
Davies and Walters measured evapo- 
ration from porous wet plates placed 
in the free air stream. The highest 
air velocity in any of the tests was 
approximately 10 m. per second, as 
shown in Table 1, which summarizes 
the conditions of these investigations. 


TABLE 1 


EXPERIMENTAL CONDITIONS UNDER WHICH DATA SHOWN IN FIGURE 1 WERE OBTAINED 


Tunnel 
Cross-section XD 
Investigator cm x cm Form Liquid cm. 
Maisel' 15.2 x 61.0 Wetted Water 18-117 
flat surface 
Powell and 54 x 54 Wetted Water 1.8-24 
Griffiths flat surface 
Wade? i2:x6 8.9 x 8.9 cm. Carbon 8.9 
Pan tetrachloride 
Trichlorethylene 
lolnene 
Benzene 
Ethyl acetate 
Acetone 
Pasquill' 19x 19 cm. Bromobenzene 19 
Pan 
Davies and Wet plate Aniline 8.9-38 
Walters? 
Millar’ 14.6 x 11.9 11.2 x 12:cm. Water Lh 
Pan 
Powell’ 54 x 54 Wetted Water 7.5-109 
evlinder 
parallel to 
flow, 4.36 cm. 
diam. 
Present data! 2.54 x 2.54 Subliming Camphor 8.8 


cylinder 7.0 

min. diam. 

parallel to 
flow 


FLAT PLATES 
MAISEL AND SHERWOOD - WATER 
POWELL AND GRIFFITHS - WATER 
WADE — 6 ORGANIC L/QU/DS 
DAVIES AND WALTERS — ANILINE 
PASQUILL — BROMOBENZENE 


ereeee==- POWELL CORRELATION -WATER 





PARALLEL CYLINDERS 
wmwwsmmmee =POWELL CORRELATION - WATER 
VW PRESENT DATA — CAMPHOR, M<0.3 


Approximate 
Range of 
Air Speeds 


cm/sec. 


100-500 


100-900 


200-1000 


60-550 


4250-59400 














OO/ MILLAR - WATER 
0.008 
Yp 
op 7-006 
Crye2 
0.004 
TURBULENT —- PRANDTL SCHL/C 
0.002 | | | 
0.00/ 
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Figure 1—Comparison of data on mass transfer from flat plates with friction data (subsonic). 
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\iaisel’s results were obtained with 
various lengths of dry approach sec- 
tion preceding the ‘wetted surface, 
and his correlation of these data in- 
troduced a correction term to allow 
for the different positions of the wet- 
ted surface in relation to the leading 
edge of the plate. The wetted surfaces 
used by Powell and Griffiths and by 
Davies and Walters were fitted with 
very short dry approach sections. In 
the case of a pan placed flush with 
the floor of a wind tunnel the dry 
“approach sections” were quite long, 
and the flow over the test surfaces 
was fully turbulent. For these reasons 
the data cannot be compared pro- 
perly in relation to a Reynolds num- 
ber defined in terms of a single char- 
acteristic length, but for simplicity 
Re, is used, with length taken as that 
of the wetted surface. 


Powell does not tabulate his data 
on evaporation from cylinders placed 
parallel to the air stream, but his data 
points fall quite close to the dashed 
line shown. Within the precision of 
the general correlation, this agrees 
with the data on flat surfaces. Because 
of this, and since both Sogin and 
Jacob@) and Chapman and Kester!) 
found little effect of surface curvature 
on heat transfer and skin friction, the 
results of the present study using 
parallel cylinders are assumed to be 
comparable with data on flat plates. 


The four points of Figure 1 at the 
highest values of Rep are data ob- 
tained in the present study at the four 
lowest air speeds employ ed. They are 
for the sublimation of camphor from 
a coating 4.6 cm. long on a parallel 
cylinder 7.0 m.m. in diameter. 


The evaporation data may be com- 
pared with the skin friction coeffici- 
ents, represented by the solid line. The 
straight section at the left (RE < 6000) 
represents C,,/2 for laminar boundary 
layers, and the slightly curved right 
portion (RE>6000) is the Prandtl- 
Schlichting relation for C,/2 in a 
turbulent boundary layer. These are 
average coefficients based on the total 
length x, from the leading edge. C; 
falls off slowly with downwind length, 
and C,/2 and jp are strictly compar- 
able only for evaporation tests in 
which the wetted surface begins at 
the leading edge of the plate. 


In view of the widely different ex- 
perimental conditions used by the 
different investigators and the omission 
of a correction for length of approach 
section, the general agreement of the 
evaporation data may be considered 
to be surprisingly good. Furthermore, 
the close check with the line for C,/2 
lends excellent support to the basic 
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Figure 2—Diagram of experimental test equipment. 
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Figure 3—Aluminum model used in Mach 2.02 wind tunnel. 


mass-momentum transfer analogy, as 
represented by the Chilton-Colburn 
representation of the effect of the 
Schmidt group. Similar data on heat 
transfer between flat plates and sub- 
sonic air streams have been shown by 
McAdams‘!2) to support the equality 
of jy and C,/2. 


Experimental Procedure 

The technique used was to immerse 
a streamlined test piece of a sublim- 
ing solid in a high-velocity air stream 
and to determine the rate of sublim- 
ation both from the weight loss and 
the change in dimensions. Figure 2 is 
a diagram of the general arrangement 
of test piece in the high-velocity 
nozzle. Figure 3 describes one of the 
three test pieces used. This is of alu- 
minum, with a cut from B to C which 
is filled with the subliming solid to 
make a smooth cylinder 10.3 m.m. 
diameter. The boundary layer trip in- 
dicated at A was employ ed to ensure 
a fully turbulent boundary layer over 
the test surface. 


The test models were prepared by 
repeatedly dipping the metal cone- 
cylinder into the molten compound 
until the whole was covered. The 
model was then secured in a lathe and 
the organic material machined until 
a smooth surface was obtained over 
the section BC of the same diameter 
as the metal ends, with no excess 
material on the metal surfaces. The 
diameter was measured accurately at 
three 120° positions at points 10 mm. 
apart along the section BC. The 
whole model was then weighed and 
positioned in the air stream for periods 
of 9 to 91 minutes. It was then re- 
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moved, weighed, and the triple mea- 
surements of diameter at 10 mm. 
intervals along the length repeated. 


Three separate tunnels were used 
to explore the velocity range from 
Mach 0.12 to 2.02. The first, used in 
the subsonic tests, exhausted directly 
to the atmosphere. The test section 
had a small cross-section (2.5 x 2.5 
cm.) so the boundary layers on model 
and wall occupied a small but appre- 
ciable fraction of the clearance. Static 
and Pitot tube pressures were mea- 
sured, as well as surface, ambient air, 
and stagnation temperatures. The last 
employed a thermocouple placed 1 
mm. inside a small impact tube. Thy- 
mol, camphor, and naphthalene were 
sublimed in these subsonic tests. 

The second was a small variable- 
Mach number recirculating tunnel of 
a larger cross-section, operated from 
Mach 1.1 to 1.9. Camphor and naph- 
thalene were sublimed in this equip- 
ment. Static pressure was measured 
at the test section and stagnation pres- 
sure and temperature at an upstream 
low-velocity reservoir. Temperature 
of the model surface was measured 
by thermocouples. 

Similar tests at Mach 2.02 were 
made using naphthalene in a larger 
recirculating tunnel having fixed 
nozzle blocks and a test section 20 x 
20 cm. Toepler Schlieren photo- 
graphs of the model were obtained in 
each supersonic test. The Mach num- 
bers calculated from the shock wave 
angle at the tip were found to agree 
well with the values obtained from 
the pressure measurements, using the 
isentropic expansion relation 
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Figure 4—Photograph of test speci- 
men in tunnel. 


P, k-1 a 


Figure 4, taken in the second tunnel 
at Mach 1.72, is one of the Schlieren 
pictures; unfortunately the reproduc- 
tion does not show clearly the shock 
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at the tip. The shocks shown origi- 
nating at the tunnel walls were present 
only in the variable-Mach-number 
tunnel. 

Static (free-stream) temperatures T, 
of the air were calculated from M, 
by the equation 


ea pom BEE (4) 


Resulis 

The average mass transfer coeffici- 
ent in each test was calculated from 
the loss in weight AW in time 0, 

using the equation 
AW e 
kc = (5S) 

m,0A (Pw — Pro) 

where A is the surface area and m, 
the molecular weight of the test ma- 
terial. Since the maximum weight loss 
(at Mach 2.02) was 2.5 g. and the 
diameter decrease seldom more than 
1.0 mm., a mean value of A was used, 
and P,,, could be taken at zero, even 
in the recirculating tunnel. At these 
low rates of evaporation the “wet- 
bulb” cooling of the surface was quite 
negligible. The driving force, there- 
fore, was taken to be the vapor pres- 
sure of the material at the temper- 
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Figure 5—Typical results showing j, based on local values of mass transfer 
coefficient. 
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atures measured at the metal cylinder 
surface. 


Local coefficients were obtained at 
several positions along the test cylin- 
der from the measured changes in 
diameter and the known density of 
the solid. Figure 5, in which jp is 
plotted vs. Re,, shows a portion of 
the results obtained in this way. Here 
Jp involves the Jocal value of kg for 
the 1.0-cm. test length, and Re, is 
based on the distance x from the tip 
to the point of measurement. The 
solid line represents the skin friction 
curve (C;/2) corresponding to 
Prandtl-Schlicting /ocal values of the 
skin-friction coefficient. 

These data suggest two conclusions: 
(a) jp is considerably greater than 
C,/2 at the lower air speeds, and (b) 
there is a clear decrease of j) with M, 
at a fixed Reynolds number. 

Average values of kg and jp were 
obtained both by averaging the local 
coefficients and from the measure- 
ments of weight loss. In some tests, 
particularly at the highest air speeds, 
there was a tendency of the sample to 
erode noticeably at the front edge 
near point B of Figure 3. Although 
no flaking or other mechanical loss 
was observed, any such loss would 
introduce a second experimental error. 
Since both of these errors would in- 
fluence the weight loss more than the 
average of the numerous dimension- 
change measurements, the average of 
the local coefficients is used in the 
subsequent discussion of the results. 
The two averages actually agreed 
within 15% in the majority of the 
tests. 


The effect of air velocity on kg is 
better shown by Figure 6, in which 
the — dimensionless group 
involving kg is plotted vs. Rey, the 
latter being based on the total length 
from trip to downstream end of the 
test section. Since the slope is nearly 
unity the variation in x, (or the 
choice of xp or Xp) is not important. 
The only other data on evaporation 
from a parallel cylinder, those of 
Powell at low air speeds, are repre- 
sented by the dashed line, and a best 
line through the mass of points on 
Figure 1 is represented by the dotted 
line. The subsonic sublimation data 
fall along a single line, but the super- 
sonic points fall as much as 50% 
lower at the same Re». 


Above Mach 1.0 the rate of evapo- 
ration is evidently reduced by some 
effect of gas compressibility. Put an- 
other way, the traditional correlation 
of mass and heat transfer cofficients 
in terms of mass velocity is not use- 
ful—changes in gas velocity and gas 
density influence kg in different ways. 
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Ejfect of Mach Number 

[hat the rate of evaporation should 
be influenced by variations in Mach 
number at a fixed Reynolds number 
is to be expected from the analogy 
to friction, where such an effect is 
well known. Kinetic energy and en- 
thalpy are interchanged in the boun- 
dary layer, and in the case of an in- 
sulated flat plate the sum of the two 
is essentially constant from free stream 
to plate surface (the “recovery factor” 
is 0.89). As a result of the temperature 
gradient, the fluid properties vary 
across the boundary layer. The insu- 
lated plate is at a higher temperature 
than the free stream, and the friction 
is less than would be predicted on the 
basis of the free-stream properties. 

Because of its importance in high- 
speed flight, this effect of Mach num- 
ber on skin friction has been studied 
by many investigators. The results 
of the various theoretical and experi- 
mental studies are well reviewed by 
Shapiro") and by Rubesin, Maydew, 
and Varga‘!4), The nature of the in- 
dicated Mach number effect is shown 
by Figure 7, which shows curves 
based on the theories of von Karman 
(13), Frankl and Voishel(5), and Van 
Driest‘'®), The various theories pre- 
dict a decrease in skin friction of 10 
to 35% at Mach 2, not depending 
greatly on Rey. 

The experimental data on skin fric- 
tion are found !9)(14) to be in good 
agreement with the Frankl and 
Voishel theory, as corrected and ex- 
tended by Rubesin, Maydew, and 
Varga. Furthermore, Rubesin(“!7) and 
others(18) (19) (20) have shown the heat 
transfer and skin friction coefficients 
to be affected by M, in the same way. 
Thus the Frankl and Voishel curve of 
Figure 7 is reasonably well established 
as representing the effect of M, on 
both skin friction and heat transfer. 

The data of the present study are 
not sufficiently precise nor suffici- 
ently extensive to establish the corres- 
ponding relation for mass transfer. 
It is evident, however, that the effect 
of M, is similar. Figure 8 shows the 
sublimation data obtained at super- 
sonic speeds, plotted as y vs. M,, 
where y is the ratio of observed k, to 
k, for incompressible flow at the 
same RE, (i.e. the ratio of the or- 
dinate of the point to that of the solid 
line of Figure 6). The solid line as 
drawn on Figure 8 indicates a reduc- 
tion in mass transfer rate of more 
than 50% at Mach 2.0, at which the 
Frankl and Voishel theory would 
predict only a 24% reduction. 

The indicated effect of Mach num- 
ber on mass transfer appears abnor- 
mally large, and may well be exag- 
gerated by experimental uncertainties. 
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The technique of sublimation at high 
speeds is new, and needs further de- 
velopment. Not only are the individual 
tests at high Mach numbers subject 
to error, but the basis of comparison 
(the solid line of Figure 6) is by no 
means well established. The line re- 
presenting incompressible flow was 
placed by drawing a straight line 
through the subsonic points of Figure 
6, without reference to the fact that 
straight line extrapolations to high 
Reynolds numbers of the low-speed 
evaporation data of Figure 1 would 
fall somewhat lower. It appears safe 
to conclude that, within the preci- 
sion of the data, the effect of Mach 
number on mass transfer is quite 
similar to its effect on heat transfer 
and skin friction. 

It should be noted that the coordi- 
nates of Figure 6 involve physical 
properties of air evaluated at the free 
stream conditions, and that k,, based 
on the vapor pressure at the surface 
temperature, is used in place of the 
coefficient k, suggested by Equation 
|. If k, were retained, with a driving 
force based on volumetric concen- 
tration corresponding to the vapor 
pressure at the surface, the product 
k;RT, would be used in place of 
k;RT,, as actually employed in pre- 
paring the graph. This change would 
amount to multiplying the ordinate of 
Figure 6 by T,/T,, and would raise 
the supersonic data nearly to the solid 
line extending the subsonic data. With 
the further modification of replacing 
u, by the air viscosity at T, the 
correlation obtained would be quite 
good, Evaluation of both density and 
viscosity at the adiabatic wall tem- 
perature is, in fact, the basis of the 
von Karman theory represented by 
the lowest line of Figure 7. The more 
refined theories allow for variation of 
physical properties with temperature 
through the boundary layer; the use 
of the arithmetic mean of T,, and T, 
has been found empirically to cor- 
relate the supersonic data on_ skin 
friction(!3), 

Figure 9 shows the data replotted 
as in Figure 6 but with Mach number 
correction applied by introducing the 
factor y, obtained from the solid line 
of Figure 8. This serves to establish 
the supersonic data on the line for 
incompressible flow, and provides a 
reasonable correlation of the data on 
cylinders over the extreme range of 
air velocities from the 60cm. per sec. 
used by Powell to more than 50,000 
cm. per sec, in the present tests. 
Appendix 

Data from the literature on the 


vapor pressures of the test materials 
were found to be in disagreement. 
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Because of this, and because of the 
po sibility that cast and finely-divided 
cr\ talline solids might have different 
vapor pressures, an experimental study 
of the vapor pressures of camphor, 
naphthalene, and thymol, was under- 
taken, 

\ir was passed through a light alu- 
minum tube filled with broken pieces 
of the cast material, the packed depth 
being sufficient to ensure essentially 
complete s saturation of the slow air 
stream. Vapor pressures were calcu- 
lated from the measured weight loss 
and air volume passed. Values of 
vapor pressure obtained in this way 
are well represented by the following 
equations: 


2855 cm 
login py = — T + 8.926 (12 to 45°C.) 
Naphthalene: 

3765 23 wes 
logiopy = — T + 11.55 (0 to 38°C.) 
Thymol: 

4780 bs 
logiopy = — 1 + 14.24 (10 to 50°C.) 


In these expressions T is in °K. and 
p, is in m.m. Hg. 


Nomenclature 

A = area of evaporating surface, sq. cm. 

C; = local skin-friction coefficient 

Cy = average skin-friction coefficient 
over length x1 

Cp = specific heat at constant pressure, 


cal/(g) (°C.) 


D = molecular diffusivity, sq. sm/sec. 

E = eddy diffusivity, sq. cm/sec. 

h = heat transfer coefficient, g cal 
(sec. ) (se. em.) (°C.) 

jp = dimensionless ratio involving kg, 


defined by Equation 2. 


jn 


ma 


P. 
vw 


XD 
XT 





dimensionless ratio involving h, 

defined by Equation 2. 

ratio of specific heat at constant 

pressure to specific heat at con- 

stant volume (taken as 1.40 for air) 

= mass transfer coefficient, g moles 

(sec.) (sq. cm.) (g mole/cu. cm.), 

= kgRT 

mass transfer coefficient, g moles 

(sec.) (sq. cm.) (atm. ) 

molecular weight of gas stream 

(= 29) 

molecular 

material 

= Mach number of free stream 

= vapor pressure, m.m. Hg. 

= static pressure of free stream, atm. 

= stagnation pressure, atm. 

= vapor pressure at surface tempera- 
ture, atm. 

= partial pressure of evaporating 
substance in free stream, atm. 

cm.) (atm)/ (g 


weight of subliming 


ll 


= gas constant, (cu. 
mole) (°K). 

= absolute temperature, °K. 

= temperature of free stream, °K. 

stagnation temperature, °K. 

= surface temperature, °K. 

free stream velocity, cm./sec. 

distance downstream from leading 

edge of plate or tip of ‘cone- 

cylinder, cm. 

length of evaporating section, cm. 

= distance from leading edge to 
downstream boundary of evapo- 
rating section, cm. 


ll 


Dimensionless groups: 


Pr = Prandtl number 

Rep Reynolds number = VopPoxp/Mo 
Rex Reynolds number = Vopox/ Mo 
Rer Reynolds number = Vopoxr/ Mo 
Greek: 

€ = eddy viscosity, sq. cm./sec. 

0 = time, sec. 

Le = viscosity, g/(sec.) (cm.); Mo at To 
p = density, g/cu.cm.; po at To 

y = mass transfer rate with supersonic 
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flow as a fraction of the rate in 
subsonic (incompressible) flow at 
the same value of Rer. 
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The Upwards Vertical Flow of 
Air-Water Mixtures 


I. Effect of Air and Water-Rates on 
Flow Pattern, Holdup and Pressure Drop 


G. W. GOVIER?, B. A. RADFORD*, and J. S. C. DUNN* 


The elementary thermodynamics and the flow pattern terminology 
of the upwards vertical flow of gas-liquid mixtures is reviewed and new 
terminology is proposed. 

Data on the upwards vertical flow of air-water mixtures in a smooth 
bore 1.025 in. I.D. tube are presented for a range of nine water rates 
from 0.00040 to 0.0421 cu. ft. per see. and for air-water volume ratios 
under flow conditions from 0 to 348. The properties of the gas and 
liquid phases were held constant through the maintenance of constant 
average flowing pressure and temperature at 36.0 psia. and 70°F. 

A correlation of the data based upon the thermodynamic analysis 
is presented. The correlation enables the prediction of flow pattern, 
pressure drop and holdup ratio for gas and liquid rates within the range 
tested but is restricted as to tube diameter and liquid and gas phase 
properties. 


Introduction amounts of the phases, their 
dispersion one in the other, and 
their individual motions are all 


important. 


HE general problem of the vertical 
flow of gas-liquid mixtures has 
received attention by a number of 
workers from the early days of the 
invention of the airlift pump to the 
present time. The application of gas 


(2) The two phase mixture actu- 
ally in motion in a vertical tube 
is never of the same overall 
composition as that admitted or 


lift principles to oil well production 
in the late 1920’s stimulated further 
interest in two phase vertical flow and 
resulted in a number of mathematical 
and experimental studies. These inves- 
tigations have shown that the vertical 
flow of two phase mixtures differs 
materially from that of single phase 
systems in two important respects: 
(1) The flow pattern cannot be 
described simply (as for single 
phase systems) as laminar, tran- 
sitional or turbulent. In the 
complex flow pattern of two 
phase systems the _ relative 


released from the tube even 
under conditions of steady 
state. This difference in relative 
quantities of the two phases 
has been referred to under the 
headings of “slip” and “holdup”. 
In general the fluid of lower 
density tends to slip past that 
of higher density so that, in the 
case of a gas-liquid system, the 
gas-liquid ratio in the flow tube 
is reduced over that of the 
entering or leaving mixtures. 


The complexity of the flow pattern 
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and the existence of holdup make it 
obvious that the relationships between 
pressure drop and the pertinent vari- 
ables will also be much more complex 
than is the case in the flow of single 
phase fluids. Moreover the differences 
in terminology used, the many vari- 
ables involved, and the lack of or 
differences in the control of certain 
of the variables has made the inter- 
pretation and correlation of experi- 
mental data difficult. While consider- 
able progress has been made in the 
overall and qualitative understanding 
of the problem much remains to be 
done from a quantitative viewpoint. 
The work to be reported in the 
series of papers of which this is the 
first was undertaken with the purposes 
of 

(1) clarifying the terminology and 
descriptions of the flow pat- 
terns encountered, 

(2) obtaining data, under controlled 
conditions, on the influence of 
the pertinent variables on flow 
pattern, holdup and_ pressure 
drop, 

and 

(3) obtaining a general correlation 
which would enable the pre- 
diction of flow pattern, holdup 
and pressure drop for any gas- 
liquid system in upwards ver- 
tical flow. 

This paper includes a thermodynamic 
analysis of the problem and deals 
specifically with the upwards vertical 
flow of air-water mixtures. Experi- 
mental data and certain correlations 
are presented on the effect of air- 
water volume ratio, at various water 
rates, on the flow pattern, the holdup, 
and the pressure drop. 


Variables Involved 
The primary variables involved in 
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the upwards vertical flow of gas-liquid 
mixtures are the following: 


Liquid Properties 
Density, p;, or specific volume, v;, 
Viscosity, Uy, 





Surface tension (to gas), 6 


Gas Properties 
Density, pcg, or specific volume, vg 
Viscosity, Uc 

Geometrical Properties 


Tube diameter, D 
d 


D 


CONSTANT WATER RATE 





Tube roughness, 


Tube height, X h FLOW PATTERN 


BUBBLE _ ie pe a 


f 
7 Si U64—=———|— FROTH a RIPPLE b= |FILM— == 
Flow Rates 


Mass flow of Liquid phase, L 
Mass flow of Gas phase, G 


Certain of these basic properties are 


Ys 0.135 ft. per sec. 





dependent upon other variables. For 
example the gas phase density will 
depend upon the nature of the gas. 
the extent of vaporization into it of 
the liquid, the pressure and the tem- 
perature. Similarily the liquid density, 
viscosity and surface tension will de- 0 100 200 300 400 


UNIT PRESSURE DROP, ySE, ft. of water per ft 





pend upon the nature of the liquid, AIR-WATER VOLUME RATIO, R,, cuft. per cu ft 


the temperature, and the solubility in 


it of the gas— —this latter depending Figure 1—Typical pressure drop relationship at low liquid rate. 


again on pressure and temperature. 
Despite the complex interrelationship 
of properties with pressure, temper- 
ature and solubility it seems clear 


PRESSURE 
REFERENCE REGIME 


| | | 
GOSLINE (1) F}—-GAS——-++~- GAS PISTON — ——_-}>——_ LIQUD 
1936 DISPERSED | DISPERSED 


CROMER and r \? 

HUNTINGTON (8) BUBBLE} PISTON and SLUG+~ 
1940 

et al la” |. wo x —+IQUID . TURBULENT — GAS TURBULENT 
1944 viscous _| Vance! 








BERGELIN ANNULAR —— 
1949 





RADFORD SLUGGING MIXED 
1949 | FROTHY 





CALVERT and 


WILLIAMS (5) AERATED += —— PISTON ——=1=— CHURN —=1—— WAVE ANNULAR 
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GALEGAR 
et al (4) 
1956 
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Figure 2—Interrelationship of flow pattern terminology. 
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that the effects of pressure, temper- 
ature and solubility need not be con- 
sidered as separate from those of the 
basic properties listed above. One 
may, therefore, expect that the pres- 
sure drop per unit height, the flow 
pattern, and the slip or holdup would 
all be functions of the primary vari- 
ables. i.e. 


dP Z 
dX 


6 
¢d’ (w. ML, O, Vc, Mc, L, 4, D, >) (1) 


Flow Pattern = 
) : (2) 
'D 2 


6 3 
D (3) 


These relationships point out a diffi- 
culty which has not always been ap- 
preciated by other workers. It is this: 
to the extent that the specific volume 
of the gas phase (in particular) in- 
fluences flow pattern, holdup and 


t 


” { vi, ML, oO, Vc, Me, L, G, I 


Holdup = 


g”’ (w. ML, 9, Vc, Mc, L, G, D, 


60 











} 
i 
i 
' 





‘~~ anenein atte ennasettienssonennaninnnnancaniisinoseeuencateegisinern 


2s ? 232 233 237 239 242 
oven —— i iv on y a 
SLUG 3 TART OF FROTH FROTH RIPR, £ rum 
3-——Flow patterns and pressure drop regimes. 


pressure drop, this variable should be 
controlled in any experimental study 
of the effect of the other variables. 
The most practical approach to this 
seems to be through the conduct of 
tests at constant average temperature 
and at constant average pressure. An- 
other point clarified by the above 
equations is that, in principle, one 
should be able to correlate the pres- 
sure drop with the pertinent variables 
indicated without separately includ- 
ing flow pattern or holdup. Certainly, 
however, a pressure drop correlation 
would be expected to reflect the 
basically different flow patterns. 
While it seems reasonable the sur- 
face tension and liquid and gas vis- 
cosities would be variables, work of 
Gosline™) and Poettmann and Car- 
penter‘?) in particular suggests that 
the effect of surface tension is small 
and that, in regions of turbulent mo- 
tion, gas phase viscosity may not be 
important, 
Thermodynamic Analysis 


Consider a differential height, dX, 
of a smooth vertical tube of diameter, 


D, to which is admitted a liquid and a 
gas phase at rates L and G. Assume 
the phases are saturated with respect 
to one another upon admittance. Thus 
for the time interval dT the well 
known mechanical energy balance 
may be written for each phase 


Vi. dV 
Lat (map +ax-+ ES ‘4a,)=0 
Ke 
(4) 


V'c dV'c 
“+ dF = QO 
Le 


(5) 


GdT (woar+ dX+ 


where 


L- = rate of supply of liquid, lb/hr. 
3 rate of supply of gas, lb/hr. 


ll 


v_ = specific volume of liquid, ft.*/Ib. 

VG = specific volume of gas, ft.3/Ib. 

dX = differential height, ft. 

P = pressure, lbs/ft.? absolute 

V’. = actual linear velocity of liquid, 
ft./sec. 

V’c = actual linear velocity of gas, ft. 
sec. 

dF, = irreversibility effect attributable to 


liquid in height dX, ft. liquid 

dFg = irreversibility effect attributable to 
gas in height dX, ft. gas 

dT = differential element of time, hrs. 
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associated with interfacial surface.) CONTROLLER 7) 
Adding these relations yields 
Vii dV'L 
vidP + dX = + dF. + 
Le Z| 
PLUG VALVE | 
G ap + Sax 4 GV eds | | 
Co ae L g. DOWNSTREAM | 
C PRESSURE STATION | 
x 
, dFo = 0 (6) | AIR- WATER 
? | SEPARATOR 
If the kinetic energy effects are | 
neglected, or considered to be in- | 
cluded within the terms dF,, and dFg, WATER ROTAMETERS | 
then - MID—POINT DRAIN 
; ae PRESSURE STATION 
G . G 2B j (TO PRESSURE ee 
dP { vi + 1 Ye i SE a I 7 a CONTROLLER ) 
4 “ 
(dF) =O (7) | 1" AIR SUPPLY 
where | 75 PS.1.6. 
G | VALVE LINK Fi 
(dF). = dF. + 7 dFo (8) | 
< . FILTER | 
Rearranging, one obtains oa | 
4 Ee 
dP fish Rea 1 dF \"" WATER SUPPLY UPSTREAM AIR ORIFICE 
VL ™ 90 PS.I.G | PRESSURE METERING 
dX 1+ Rv 1+Rv\dX/, STATION STATION 
where PLUG VALVE 
oe try x. & 
Rm = gas-liquid mass ratio, I 
‘ui ., GYG DRAIN 
Ry = gas-liquid volume ratio, I (9) 
VL 


Figure 4—Schematic flow diagram of equipment. 


This equation indicates that the 
total unit pressure drop, measured in 1.0 
terms of feet of the flowing liquid 


aP . 
is 
L r? 
dX 
equal to the sum of two components 


1+ Ru 


1 + Ry 


phase per foot of height, —v 
0.2 
, the hydrostatic head com- 


ponent, ft. of liquid per ft. of height, 
and 


, ft of water per ft 























1 dF oe ate alx 06 x | « 
- the irreversibility BG x 
1+Ry \adx/, ' | | ; 
i s 
sec CONSTANT WATER RATE 
component, ft. of liquid per ft. of “ V, © 0.0698 tt. per sec. = 
height. Oo ps 
Equation (9) may be applied to a e ‘ae 
system of finite height only if the ' a 
fractional variations in R,, and Ry are Ww = 
small enough to justify the use of = > 
average values based upon an average us 
system pressure and temperature. x 
Under these conditions 0.2 
J ad 
AP 1+ Ru 4 1 AF S 
-. i re > 
a ith THRE CO 
(10) 
is ) 1 oS oO , ‘ p CX- 
(It is of interest to note that the ex ‘ - — — <a a 
, AP. . 
pression — y, 3 is exactly equivalent AIR-WATER VOLUME RATIO, R,, cu. ft. per cu. ft 
‘ eis Figure 5—Pressure drop and holdup relationship, V,—0.0698 ft/sec. 
to the “submergence ratio” referred : 
to in works on air lift pumps.) 
° ee ee > . . 2f, V2 
In the limiting condition when eae SP ee AFL (11) AFL _ 2fiVi (12) 
Ry=Ry=O, or a single phase liquid AX AX AX g-D 
is flowing, the equation reduces to where and f, is the common friction factor 
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DVL 6 Scat 
dependent upon and 5 AF 2’, V3, .; 
MLVL Z = (13) 
AX}, g- D 


In the general case a two phase 
friction factor f’;, may be defined 
through the relation 


where V,, is the superficial liquid 
phase velocity based upon the total 


tube cross section. This two phase 
friction factor might be expected to 


DVL ‘ 
depend upon ,. and either the 


VLEL 
mass or the volume ratio of the two 
phases. 

An alternate rearrangement of equa- 
tions (4) and (5) yields 


L L 
dP € vie + “) + dX € + ') Fr 


(dF); =O (14) 
where 


L 
(dF )c = — dF, + dFc 
G 
Under conditions where average 
values of R,, and Ry may be used 
this leads to 


1 
1+ 
AP Ru + 1 AF 
-— Vo = 
AX 1 f VAR, 
1+ 1 + : 
Ry Ry 


(16) 
relating the total unit pressure drop 
measured in terms of feet of the flow- 
ing gas phase per foot of height, 

P 

ae VG = 
x 
static and irreversibility components. 


, to corresponding by hydro- 


In limiting condition of the flow of 


1 
a single gas phase = = O and 
eee a , 
the equation reduces to 
AP 4 AF« (17) 
= Vo = / 
AX AX 
where 
AFc 2fcVc? 
Ne (18) 
AX g- D 


and f,, is the common friction factor 
DVe 6 
dependent upon and 
UGVG I 
Again, in the general case, a two 
phase friction factor f’, might be de- 
fined by the relation 


AF 2f’cVc? 19) 
Ans... &D ( 


where V, is the superficial gas velo- 
city based upon the total tube cross 
section. The choice between equa- 
tions (10) and (13) in terms of ft. of 
liquid, and equations (16) and (19) 
in terms of ft. of gas would appear 
arbitrary as the two are mathematic- 
ally equivalent. 

These basic relationships suggest 
that pressure drop data for the ver- 
tical flow of two phase systems might 
best be understood 

(1) by separating the unit irrever- 

sibility component from the 
total pressure drop, 

(2) by calculating a two phase 
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friction factor (either f’;, or f’, 
from equation (13) or (19)) 
and 
(3) by seeking a correlation of the 
friction factor with a super- 
ficial Reynolds number, a rela- 
tive roughness term and a mass 
or volume ratio of the phases. 
The fact that even a very small 
amount of liquid in a gas stream 
greatly increases pressure drop indi- 
cates further that the presence of the 
liquid phase may well outshadow the 
effect of roughness‘ 3) and also that at 
high velocities (fully developed tur- 
bulence) viscosity effects may be very 
small. In this work, where tube dia- 
meter and gas and liquid properties 
are held constant, the pressure drop 
data are correlated as functions of 
superficial liquid velocity (propor- 
tional to liquid Reynolds number ) 
and the volume ratio of the phases. 
Poettmann and Carpenter?) have 
also applied the mechanical energy 
balance to the vertical two phase flow 
problem but in an approximate man- 
ner. Their assumption that a two (or 
three) phase mixture could be treated 
as a pseudo homogeneous phase of 
constant equivalent specific volume 
has led to a useful practical correla- 
tion but is not too sound fundamen- 
tally. 


Hold-up and Slip 


As has been mentioned one of the 
distinguishing characteristics of two 
phase vertical flow is the develop- 
ment, within the flow section, of 
concentration of the denser phase 
greater than that in the supply mix- 
ture. The term “holdup ratio” is used 
in this paper as a quantitative measure 
of this liquid accumulation. Holdup 
ratio is here defined as the ratio of 
(a) the gas-liquid volume ratio in the 
supply mixture to (b) the gas-liquid 
volume ratio in the flow section. Both 
volume ratios are expressed in cu. ft. 
per cu. ft. at average flow pressure 
and temperature. Galegar, Stovall and 
Huntington'*) have reported quanti- 
ties equivalent to the individual ratios 
(a) and (b) and Calvert and Wil- 
liams(5) have made a mathematical 
analysis of holdup i in the film or a. 
lar flow region. Other workers(!: +5.7 
have referred to the “slip velocity” Ne 
the gas phase. Slip velocity has been 
defined as the difference in the aver- 
age lineal velocities of the gas and the 
liquid phases both measured at aver- 
age flow pressure and temperature. 

The two quantities, holdup ratio 
and slip velocity are two measures of 
the same phenomenon and are simply 
interrelated as follows: 

Let the volume fraction of the flow 
section occupied by the gas phase be 
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Figure 8—Pressure drop and holdup relationship, 


°(See Note, page 22). 
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Figure 9—Pressure drop and holdup relationship, 


°(See Note, page 70). 


defined as q. Then 1-q represents the 
volume fraction occupied by the 
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Figure 10—Pressure drop and holdup relationship, 


*(See Note, page 70). 
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Figure 11—Pressure drop and holdup relationship, V,=3.47 ft/sec. 


*(See Note, page 70). 


Let A be the cross sectional area of 


the empty flow section. Then the 
average lineal velocities are 
—" Gvc 
\ G = 
Aq 
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i, Lyi 
A A 
and the slip velocity = 
1 / Gvc Lyi 
A\ q “i 


—_ q) 
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HOLDUP RATIO, He 


Hp 


HOLDUP RATIO, 


Lv 
Vo —- Vit = (Hr — 1) = 
A (1 — q) 
V't (Hr —.1) 
or 
H Jan: \ 's _ V1 4 1 
R = Vv, 


Flow Pattern Terminology 


The flow pattern for the upwards 
vertica! flow of gas-liquid mixtures 
has been variously described by Gos- 
line“), Bergelin), Stovall and Hunt- 
ington’), Calvert and Williams), 
Cromer and Huntington *), Radford 
(9), Dunn 9), Martinelli, mR. ‘C. et al. 
(11) and others. The descriptions of 
most of these workers were based 
either upon visual and qualitative 
studies or upon “efficiency” or liquid 
“vield” considerations. Radford and 
Dunn and, later, Calvert and Williams 
correlated their flow pattern descrip- 
tions with pressure drop observations. 
This method, which has the advantage 
of greater precision of definition, is 
used here. 

A typical pressure drop vs gas- 
liquid volume ratio curve at constant 
low liquid rate and at constant aver- 
age pressure and temperature for air- 
water flow in a smooth tube of 1.025 
in. id. is shown in Figure (At 
higher liquid rates the curve changes 
its shape and shifts upwards as will be 
discussed later.) On the basis of this 
relationship it seems reasonable co 
divide the curve as indicated into four 
pressure drop regimes I, II, iff, IV 
characterized respec tively by a de- 
creasing, an increasing, a decreasing 
and an increasing unit pressure drop 
with increasing gas-liquid volume 
ratio at constant liquid rate. This basic 
division may then be subdivided, and 
given qualitative description by re- 
ference to visual observations. The 
terminology used by previous investi- 
gators has been reviewed and com- 
pared on this basis. Figure 2 shows 
the interrelationship of ‘the terminol- 
ogy as nearly as this can be deter- 
mined from the published work. 
While the importance of flow pattern 
as determined visually is recognized, 
the primary subdivision on the basis 
of pressure drop regime is considered 
more precise and to offer greater 
possibilities as an aid to correlation. 
The flow patterns encountered in the 
four pressure drop regimes are illus- 
trated by diagrams and photographs 
in Figure wi 


Low Liquid Rates 
Regime I 
As the gas rate, at any constant 


liquid rate, ‘is increased from zero to 
low values (gas-liquid volume ratios 
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up to 1-5) the gas is dispersed as dis- 
crete bubbles which increase steadily 
im number and size. This is the bubble 
flow pattern and it extends only a 
short way into Regime I. W ith a 
further increase in gas flow some of 
tlie bubbles coalesce to form larger 
lnullet shaped slugs. These slugs grow 
in length and diameter and their up- 
ward velocity increases as the gas 
rate is further increased. This slug 
tlow pattern is characterized by alter- 
nating slugs of gas which are sur- 
rounded by a thin liquid annulus and 
interspersed with regions of bubble 
tiow. As a slug passes any given tube 
section some of the liquid falls down- 
ward into the bubble filled region 
beneath it. Slug flow continues to the 
end of Regime I and well into Regime 
IL. 

Regime II 

As the gas-liquid volume ratio is 
still further increased the tendency 
of the liquid phase to fall downwards 
beneath the slugs decreases and finally 
stops. At this latter point the slugs 
seem to merge and degenerate into a 
froth. Froth flow would seem to ap- 
pear first at a point roughly mid-way 
in Regime II. This pattern is highly 
turbulent with coarse agitation of the 
liquid phase being its most predom- 
inant characteristic. It extends through 
the second half of Regime II into 


Regime III. 


Regime III 

The flow pattern in Regime III is 
primarily froth. As the gas-liquid 
volume ratio is increased the amount 
of relative motion between phases 
seems to decrease. This phenomenon 
continues with increasing gas rate 
until a region of ripple flow is entered. 
Ripple flow is characterized by an 
upwards moving wavy or ripply layer 
of liquid on the tube wall. The zone 
of ripple flow is rather narrow and 
actually forms a transition region be- 
tween Regimes Ill and 1V. 


Regime IV 

The ripple flow pattern continues 
a short way into Regime IV and then 
develops into film flow. In this flow 
pattern the liquid flows as an annular 
film and the gas flows at a higher 
velocity through the inner core. This 
flow pattern results from the smooth- 
ing of the waves or ripples of ripple 
flow by the attrition of the higher gas 
rates. Calvert and Williams‘) have 
analyzed the film flow pattern from 
the viewpoint of the actual shearing 
forces involved. Eventually, at some 
sufficiently high gas rate, a sufficient 
number of water droplets will be 
captured from the liquid interface 
to remove the liquid layer and to 
result in a region of mist flow. The 
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Figure 14—Minima and maximum pressure drop loci (vs R,,) 
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Figure 15—Pressure drop components. 


occurrence of this zone is obvious 
but has not been directly observed 
in this work. 


High Liquid Rates 


The above description of flow pat- 
terns does not apply in its entirety 
to higher water rates. It will be shown 
that as the water rate is increased 
Regimes Il and III become progres- 
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sively narrower and finally disappear. 
At these higher water rates the flow 
pattern of the liquid phase seems 
much the same as at the lower rates 
but the gas phase pattern changes 
markedly. At the close of Regime I 
the gas phase, as before, is character- 
ized by an agitated pattern which 
extends the entire length of the tube 
with fast moving slugs travelling in- 


termittently through the system. As 
the gas rate increases these slugs be- 
come faster moving and more fre- 
quent in occurrence. No froth flow 
pattern is observed. The amount of 
turbulence present in the liquid phase 
decreases with increasing gas rates 
and develops directly into a pattern 
similar to ripple flow. Ripple flow 
gives way to film flow as at the lower 
water ‘rates. 


Experimental Equipment 

The. experimental equipment _ is 
illustrated schematically in Figure 4. 
It consisted of an air supply and 
metering system, a water supply and 
metering system, the two-phase test 
section and a manometer system for 
measuring pressure drop. 

Air was supplied at approximately 
100 psi. from a central compressor, 
passed through a reducing valve to 
maintain a nearly constant down- 
stream pressure of about 75 psi. and 
then through an_ orifice metering 
station to the base of the test equip- 
ment. Control was effected manually 
by a suitable valve in the main one 
inch line and by a needle valve in a 
quarter inch by-pass line. Water was 
supplied through a 70 gallon per 
minute capacity centrifugal pump 
having a maximum head of about 110 
psi. The water was admitted to the 
base of the equipment through one of 
four rotameters connected in parallel. 
These rotameters had maximum capa- 
cities of 0.000272, 0.00247, 0.0211 and 
0.0219 cu. ft. per sec. thus permitting 
relatively accurate metering of the 
water over the full range from 0.0004 
to 0.0421 cu. ft. per sec. 

Preliminary tests were conducted 
with different designs of mixing sec- 
tions at the base of the two-phase flow 
column. It was found however that 
the design of the mixing section af- 
fected the flow pattern only for a 
very short distance above it and that 
with an adequate calming section no 
special mixing device was required. 
Accordingly the air and water were 
combined in a simple tee at the base 
of the test section. 

The test section itself consisted of 
a transparent section of cellulose 
acetate-butyrate plastic tubing ap- 
proximately 30 ft. long and of one in. 
nominal diameter. The lower end of 
the plastic tube was connected to the 
mixing tee with a two ft. section of 
one in. galvanized iron pipe contain- 
ing a full flow plug valve. The upper 
end of the plastic pipe was similarly 
connected to the discharge piping of 
the system through another full flow 
plug valve. These plug valves ope- 
rated with a 90° turn of the handle. 
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The handles of the two were linked 
together with a length of one in. 
aluminum angle iron. This permitted 
the two plug valves to be opened or 
closed simultaneously. Three pressure 
taps were made in the plastic tube. 
The lower or upstream tap was ap- 
proximately eight ft. above the mixing 
chamber and the upper or down- 
stream tap 22.88 ft. above it. The 
third pressure tap was located mid- 
way between the upstream and down- 
stream ones. All pressure taps were 
made in blocks of lucite, machined to 
fit around the tube and sealed to it 
with chloroform. The block and tube 
were then bored with a 3/16 hole 
which was carefully reamed and the 
block was counter-bored and tapped 
for 1/8 in. standard pipe fittings. 
Connected to the discharge end of the 
test section were an air-water separ- 
ator which permitted the water to be 
discharged to sewer and a pneumatic 
control valve installed in the air dis- 
charge line from the separator. 

Two 30 in. Meriam manometers, 
one filled with water, the other with 
mercury, were provided for measur- 
ing the pressure drop across the air 
metering orifices. A 60 in. mercury 
manometer and a 0 to 50 psig, re- 
cording Brown pressure controller 
were connected in parallel to the mid- 
point pressure tap of the test section. 
This connection was made through a 
small plastic separator directly con- 
nected to the mid-point pressure tap. 
This separator enabled the operator 
to keep the pressure lead line to the 
manometer and the Brown instrument 
clear of water by periodic draining 
from its base. The pressure controller 
actuated the pneumatic valve in the 
air discharge line from the main 
separator at the outlet end of the 
test section. Two 60 in. Meriam 
manometers were provided for the 
measurement of the pressure drop 
between the upstream and the down- 
stream pressure taps. One of these 
manometers was filled with water, 
the other with mercury. They were 
connected in parallel through plastic 
separators directly connected to the 
pressure taps. 

Thermometers were provided for 
measuring the inlet temperatures of 
the water and the air and the outlet 
temperature of the two-phase mix- 
ture. The orifice plates used for mea- 
suring the air rates, were each cali- 
brated with the aid of a critical flow 
prover. The calibrations were checked 
one against another for consistency 
between plates. The four rotameters 
used for measuring water flow were 
calibrated by a direct weighing tech- 
nique. 
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Experimental Procedure 


The tests reported here were all 
carried out at the same mid-point test 
section pressure of 36 psia. This pres- 
sure was controlled to within 0.1 or 
0.2 psia. by the Brown pressure con- 
troller. Extensive tests were made at 
each of the nine water rates. The test 
procedure was essentially the same in 
each case. 


A test was started by opening the 
two plug valves on the test section, 
closing all pressure tap connections 
and introducing water at a desired 
rate through the appropriate rota- 
meter to the base of the test section. 
With this rate set a small flow of air 
was started through an appropriate 
sized orifice to the base of the column. 
The mid-point pressure tap was then 
opened. The two-phase mixture dis- 
charging from the top of the test 
section was separated in the separator 


(1) (2) (3) (4) 


Test Water | Saturated Mid Point 


with pressure control on the entire 
system being maintained by the Brown 
instrument and the pneumatic valve 
handling the air discharge from the 
separator. With steady state condi- 
tions reached the upstream and down- 
stream pressure taps were opened, the 
small separators flushed to give an 
indication of the pressure drop 

one of the two 60 in. manometers. 
After conditions appeared completely 
stabilized, readings were taken of the 
water flow rate, the air flow rate and 
the pressure drop. A check reading 
of the mid-point pressure was taken 
on the Meriam manometer in parallel 
with the Brown instrument. Follow- 
ing this the upstream and downstream 
pressure taps were closed and the 
connecting link joining the two plug 
valves was operated rapidly to close 
both valves simultaneously and trap 
the flowing mixture in the test sec- 


TABLE 1 (PARTIAL) 


EXPERIMENTAL DaTa 


Tube Diameter 1.025 Inches 1.D. 
Test Section Length 22.88 Ft. 
Test Section Mid Point Pressure 36.0 Psia. 
Test Section Mid Point Temperature as Indicated. 


Specific Volume of Liquid Phase = .01603 Ft®/Lb. 


Specific Volume of Saturated Gas Phase as Indicated. 


(5) (6) (7) (8) 


Discharge Air 


Water Ratio** 


Test 
Section 


Specific 


Volume Holdup 


Number Rate Air Rate Temperature Air Water Saturated | Ratio*** 
Ft3/Sec. | Ft8/Sec.* °F (a) (b) Ratio Gas Phase 
Ft?/Ft®, | Lb/Lb. | Ft8/Ft®. | Ft®/Lb. 

- 00040 0. 70 0 0 0 
395D .00040 .00660 70 16.5 0483 eR 5.47 9.54 
396D 00040 0149 72 37.3 109 $3.11 5.49 12.0 
397D .00040 .0283 72 70.7 . 206 3.30 5.49 21.4 
398D 00040 0397 72 99.2 290 3.45 5.49 28.8 
399D 00040 0331 72 82.7 241 3.41 5.49 24.3 
400D .00040 0587 72 147 429 3.81 5.49 38.6 
401D 00040 0765 72 191 557 4.74 5.49 40.3 
402D .00040 139 71 348 1.02 9.47 5.48 36.7 
405D 00040 0652 70 163 477 3.86 5.47 42.2 
406D 00040 .113 70 283 829 7.10 5.47 39.9 

- 0050 0 70 0 0 0 - 
380D 0050 0620 69 12.4 .0365 2.32 5.45 5.35 
381D .0050 117 69 23.4 .0688 3.65 5.45 6.41 
382D .0050 .158 70 31.6 0926 4.34 5.47 7.28 
383D .0050 137 75 27.4 0795 4.13 5.53 6.63 
386D .0050 0303 73 6.06 0176 2.47 | 5.51 2.45 
387D 0050 | 0568 74 11.4 0332 30 | CUCU SR 5.04 
388D 0050 | .0775 72 e5..5 0451 2.81 5.51 5.52 
389D .0050 .0157 73 3.2 .00931 1.83 5.51 1.78 
390D 0050 0285 72 5.8 0169 2.28 $.5i 2.54 
391D .0050 0383 72 7.6 .0221 2.74 5.51 2.77 
392D .0050 00788 70 1.6 .00469 0.93 5.47 2.72 

- .0421 0 70 0. 0. 0. 
308S 0421 .0720 69 1.71 .00502 1.06 5.45 1.64 
309S .0421 .0360 67 .855 .00252 .628 5.44 1.36 
310S .0421 0598 67 1.42 00419 .879 5.44 1.61 
311S .0421 .0791 67 1.99 .00587 1.16 5.44 1.72 
312S .0421 .0194 72 461 .00136 .318 5.49 1.45 


* At test section temperature and average pressure. 


** Column 3 divided by column 2. 
**Column 5 (a) divided by column 6. 
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tion. Some time was allowed for the 
mixture to settle out and then the 


relative volumes of water and air 
confined between the two plug valves 
were observed with the aid of a cali- 
brated scale. Following this, condi- 
tions were adjusted to successively 
higher air rates at the constant water 
rate. Careful visual observations were 
made of the flow pattern for each 
test and several series of photographs 
were taken. 


Experimental Results 


The basic experimental data are 
tabulated in Table 1 (Partial) for the 
lowest, the highest and one of the 
intermediate water rates. (The com- 
plete tabular data including those of 
the other six water rates is given in 
Table 1 of material deposited as 
Document No. 5204 with the Ameri- 
can Documentation Institute Auxiliary 


(9) (10) (11) (12) 
Specific Volume 
Pressure Ratio Pressure Flow 
Drop V./Vg at Drop Pattern 
‘t. HoO/Ft. Paes Fay. Regime 
1.00011 I Single Phase 
0.293 29.3 x 1074 I Slug 
0.201 29.2 x 1074 I Slug 
0.195 29.2 x 1074 II Slug 
0.194 29.2 x10 Il Froth 
0.180 29.2 x 1074 II Slug 
0.198 29.2 x 1074 III Froth 
0.173 29.2 x 1074 Ill Froth 
0.108 29.25 x 1074 Ill Froth 
0.201 29.3 x107¢ Ill Froth 
0.138 29.3 x 1074 III Froth 
1.0047 - I Single Phase 
312 29.4 x107* Il Froth 
288 29.4 x 1074 III Froth 
293 29.3 x 1074 IV Ripple 
295 29.0 x 1074 Ill Ripple 
. 256 29.1 x 1074 I | Slug 
.330 29.1 x 1074 II Froth 
321 29.1 x 107% III Froth 
346 29.1 x 1074 I Slug 
272 29.1 x 1074 I Slug 
252 29.1 x 1074 I Slug 
.498 29.3 x107* I Slug 
1.19 - I Single Phase 
1.04 29.4 x 10™* IV Annular 
1.00 29.5 x 1074 I Slug 
1.02 29.5 x 1074 IV Annular 
1.05 29.5 x 1074 IV Annular 
1.05 29.2 x 1074 I Slug 
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Publications Project, Photo Duplica- 
ton Service, Library of Congress, 
Washington 25, D.C.).* 


Referring to Table 1, Column 2 
gives the water rate in cu. ft. per 
sec. Of liquid saturated gas-phase at 
the mid-point temperature given by 
Column 4 and the constant mid-point 
pressure of 36.0 psia. These rates were 
calculated from the actual measured 
air rates to correct for the saturation 
of the air with water which took 
place in the calming section. Column 
5 presents the ratio of the air to water 
rates, both on a volume and on a mass 
basis. Column 6 gives the ratio of air 
to water trapped in the test section 
at the completion of each test. Column 

gives the specific volume of the 
saturated gas-phase at mid-point test 
section conditions. Column 8 presents 
the holdup ratio, being the ratio of 
the discharge air-water volume ratio 
to the test section air-water volume 
ratio i.e. of Column 5 to Column 6. 
Column 9 is the total pressure drop 
between the upstream and down- 
stream pressure measuring stations, 
expressed in terms of ft. of water per 
ft. of height. Column 10 presents the 
ratio of the specific volumes of the 
liquid-phase to the saturated gas- 
phase at average test section pressure 
and temperature. Column 11 identi- 
fies the pressure drop regime and 
Column 12 the flow pattern. 


The basic data are summarized for 
each of the constant water rates in 
Figures 5 through 12. These figures 
show by separate curves (a) the re- 
lationship between unit pressure drop 
and the air-water volume ratio at test 
conditions and (b) the relationship 
between holdup ratio and air-water 
volume ratio at test conditions. The 
different constant water rates for each 
of the figures is represented by the 
superficial water velocity, V,. The 
existence of the four pressure drop 
regimes and the two minima points 
in the pressure drop curves is clearly 
indicated in Figures 5 through 8 
covering superficial water velocities 
from 0.0698 to 2.06 ft. per sec. and 
at least perceptible at the rate of 
Figure 9. The data for the higher 
water rates shown in Figures 10 to 12 
indicate a mergence of the maximum 
point in the pressure drop curve with 
the second minimum point. This 
mergence is reflected by the inflec- 
tion point of Figure 10 and results in 
the disappearance of pressure drop 
regimes Il and III at a superficial 
water velocity of about 3.5 ft. per sec. 


°A copy may be secured by citing the Docu- 
ment number and by remiting $2.50 for photo- 
prints or $1.75 for 35 mm. microfilm. Advance 
payment is required. Make cheques or money 
wrders payable to: Chief, Photoduplication Ser- 
vice, Library of Congress. 
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Figure 18—Holdup ratio correlation. 


The interrelationship between the 
holdup curves and the pressure drop 
curves is to be noted. In each instance 
the holdup curve exhibits a region of 
inflection about in the middle of 
pressure drop regime I and a very 
abrupt change of slope at about the 
end of pressure drop regime II. It is 
also to be noted that the holdup curve 
is nearly linear from about the middle 
of regime I to the end of regime II. 
At low water rates the holdup curve 
has a slight negative slope through 
pressure drop regime III. This slope 
gradually increases until, at a super- 
ficial lineal velocity of near 0.6 feet 
per second, the slope becomes zero. 
With increasing water rates the slope 
of the holdup curves in regime III 
increases further in a positive sense. 

The loci of the two minima and 
the maximum points in the pressure 
drop curves of Figures 5-12 are plot- 
ted in Figures 13 and 14. Figure 13 
shows the loci as functions of the unit 
pressure drop and the superficial 
water velocity. The convergence of 
the loci corresponding with the maxi- 
mum and the second minimum at a 
superficial water velocity of about 
3.5 ft. per sec. is again ev ident. Figure 
14 which usw the loci as functions 
of superficial water velocity and air- 
water volume ratio (under flow con- 
ditions) is perhaps more revealing and 
certainly more useful from the view- 
point of delineating the various pres- 
sure drop regimes and flow patterns. 
The locus of the first pressure drop 
minimum, which is the boundary be- 
tween regimes I and II, is seen to be 
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a straight line on logarithmic paper 
of slope equal to -1.17. The locus of 
the maximum point is more or less 
parallel to this and displaced toward 
higher air-water volume ratios. The 
locus of the second minimum point is 
again not too far from being parallel 
to that of the first minimum and the 
maximum points but is displaced 
further toward higher  air-water 
volume ratios. The maximum locus 
begins to approach that of the second 
minimum at superficial water veloci- 
ties above about one ft. per sec. and 
merges with it at the velocity of 
about 3.5 ft. per sec. Thus above a 
water rate of 3.5 ft. per sec. regime 
III does not exist and regime II be- 
comes continuous with regime IV. 

The loci of the points of inflection 
on the holdup curves and of the 
points of abrupt change of slope are 
seen also to be roughly parallel to the 
pressure drop minima and maximum 
loci. The approximate boundaries be- 
tween bubble and slug flow and be- 
tween slug and froth flow, froth and 
wave entrainment flow and wave en- 
trainment and annular flow are shown 
as dashed lines. These lines have been 
located on the basis of visual obser- 
vations during the tests. 


Interpretation of Results 


A considerable amount of time has 
been devoted to an interpretation of 
the pressure drop results. In the light 
of the fluid flow thermodynamics 
presented earlier the total pressure 
drop was, by calculation, separated 
into its two components—the hydro- 
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static head and the irreversibility 
component. The relationship between 
these two components and the total 
pressure drop, all in feet of water 
per foot, is indicated in Figure 15 for 
the superficial water velocity corres- 
ponding with Figure 9, namely 0.873 
ft. per sec. The components bare a 
similar relationship to the total at the 
other water flow rates. Both the fric- 
tion factors, f’;, and f’,, were calcu- 
lated for all of the test data and plots 
of these were made vs the air-water 
volume ratios at each of the nine 
constant water rates. Both sets of 
plots yielded surprisingly smooth 
curves. (Complete numerical values 
of the components and the two fric- 
tion factors are given in Table 2 of 
the material deposited as Documents 
No. 5204 with the American Docu- 
mentation Institute Photo Duplication 
Service, Library of Congress, Wash- 
ington 25, D.C.). 

Typical curves of f';, versus the air- 
water volume ratio at the lowest, the 
highest and an intermediate water 
velocity are given for illustrative pur- 
poses in Figure 16. Figure 17 was 
prepared from the set of these curves 
representing all nine water rates. This 
is a logarithmic plot of f’;, versus the 
superficial water velocity showing 
lines of constant air-water volume 
ratio. For reference purposes the 
conventional friction factor is also 
plotted against water velocity for 
single phase water flow. The friction 
factor lines representing the two- 
phase flow condition are seen to lie 
well above and to the right of the 


single-phase lines. At the lowest gas- 
liquid ratios the two-phase lines ap- 
pear to be approaching the single- 
phase line but at the highest water 
velocities fall below the single-phase 
line representing turbulent water flow. 
Clearly more data are required at 
lower Ry ratios to determine the 
manner in which the two-phase lines 
blend into the single-phase lines. 

The loci of the pressure drop 
minima and maximum points are also 
plotted in Figure 17, the points on 
each Ry curve being taken from 
Figure 14. The smooth curve loci 
resulting in Figure 17 attest the con- 
sistency of the data for the derived 
f’,, values are sensitive to minor errors 
in the total measured pressure drop. 


Figure 18 showing the holdup ratio 
as functions of air-water volume ratio 
and superficial water velocity was 
prepared from the holdup curves of 
Figures 5 to 12 inclusive. Again the 
minima and maximum pressure drop 
loci are superimposed. The relation- 
ship between these loci and the mini- 
ma and maxima in the family of hold- 
up curves is apparent. 


Conclusions 


The correlations of Equations 10 
and 13 and Figures 14, 17 and 18 
enable the prediction of the flow pat- 
tern, the pressure drop and the hold- 
up as functions of the air and water 
rates for the upwards vertical flow 
of air-water mixtures at an average 
pressure of 36.0 psia. in smooth tubes 
of approximately one in. diameter. 

These correlations apply only to a 


NOTE: Errata 


system in which the average pressure 
is 36.0 psia. and over which the frac- 
tional change in specific volume of 
the water saturated air is small enough 
to justify the use of an average value. 
Experimental work now nearly com- 
pleted and to be reported shortly has 
been undertaken to extend the appli- 
cation of the correlation to other 
diameters and other average pressures 
(or average gas phase densities). With 
the effect of gas phase density known, 
integration or step-wise use of the 
correlation will permit its application 
to systems wherein there is a large 
change in gas phase specific volume. 
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The holdup ratio curves in Figures 7-12 inclusive have 


erroneously been extrapolated to the origin of co-ordinates. 


These curves should have been extended to a holdup ratio 


of unity 


at zero air-water volume ratio. Obviously at 


extremely low air-water ratios the slip velocity is zero and 


from equation (20) the holdup ratio is unity. 
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Correlations of the Dense-Phase 
Fluidized State and their 


Applications 


MAX LEVA? 


The paper gives a brief review and application of the more firmly 
established correlations that pertain to the dense-phase fluidized state. 
It begins with a discussion of what is meant by quality of fluidization 
and then proceeds with a consideration of such phenomena as normal 


fluidization, slugging and channeling. 


Thereafter the paper reports on 


quantitative correlations that pertain to the prediction of the onset of 
fluidization, the state of bed expansion and some heat transfer character- 


istics of the fluidized dense state. 


7 purpose of this paper is to 


review some of the more firmly 
established developments that per- 
tain to the fluidized dense phase, dis- 
cuss their background and examine 
their applicability and limitations. 

Consideration will be given to: 

1. Quality of fluidization and ab- 
normalities encountered during 
fluidization. 

2. Prediction of the onset of fluid- 
ization. 


3. Prediction of the state of expan- 
sion of a column of granular 
solids. 

4. Heat transfer characteristics. 


Quality of Fluidization 

The earliest attempts to describe 
quality of fluidization were based pri- 
marily on visual observation. Since a 
fluidizing column can be observed 
only from the outside this was of 
doubtful value, especially since, due 
to the subjective nature of individual 
observations, no true standard could 
be established. This is especially true 
for fluidization of solids with gases, 
where virtually all shades of per- 
formance can be realized. The classi- 
fication aggregative and particulate 
fluidization suggested by Wilhelm 
and Kwauk ‘!) though also conceived 
by visual observation is somewhat 
more valuable. Thus aggregative 
fluidization was found to comprise 


essentially solid-gas operations, 
whereas solid- liquid operations yield 
normally to particulate fluidization. 
In the latter the particles behave 
more or less as distinct individual 
elements and thus a column exhibiting 
particulate fluidization is character- 
ized by a relatively uniform distribu- 
tion of void space throughout. This 
very lack of homogeneity of void 
space is the feature that defines 
(under the usual operating condi- 
tions) the aggregative nature of solid- 
gas systems. 

An excellent method of investigat- 
ing the internal condition of solid- 
gas systems has been suggested by 
Morse and Ballou (2). Rather than 
refer further to quality of fluidization 
they chose to define the term uni- 
formity of fluidization to—“describe 
the dispersion of particles and fluid 
in fluidized-solid suspensions”. As 
pointed out degree of uniformity 
depends on particle and fluid charac- 
teristics as well as on apparatus con- 
struction. Their method of analysis 
actually reported dielectric properties 
of sections of fluidized beds. Since 
these properties are functions of gas- 
solid dispersion, soundings could thus 
be made throughout the bed, per- 
mitting establishment of uniformity 
standards. Morse and Ballou pro- 
ceeded with some preliminary sys- 
tematic experiments that showed 
clearly that bed, fluid and apparatus 


1 Manuscript received March 8, 1957. 


2 Consulting Chemical Engineer, Pittsburgh, Penn. 
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characteristics influenced uniformity 
of fluidization. The question of uni- 
formity in fluidization is of course of 
interest to fluidization applications in 
the catalytic as well as purely physical 
and chemical fields. As of the 
moment large scale systematic 
measurements on uniformity are not 
yet available, but experiments are in 
progress ‘3) which should lead to a 
more complete and systematic under- 
standing of the factors that define and 
promote uniformity of fluidization. 


Solid-gas systems exhibiting good 
overall uniformity characteristics may 
be summarily recognized by their 
pressure drop fluid flow behavior. 
Thus in a solid-gas system of good 
uniformity pressure drops are stable 
without excessive (that is within an 
order of magnitude of = 10 percent) 
fluctuations, and moreover the value 
of the pressure drop is in good ac- 
cord with the theoretical 


Ap, = L (1 — €) (ps — pr) (1)* 


*Symbols are defined in the table of nomen- 
clature. 


Figure 1 is indicative of this be- 
havior. Branch a’ -a pertains to a 
fixed bed, and in accord with ex- 
perience of flow through packed sys- 
tems (in laminar flow), A p « G. ‘At 
point “a” the column of solids expands 
slightly, just sufficient to provide ad- 
ditional void space so that at a 
minimally higher flow rate the pres- 
sure drop can remain at the theo- 
retical value, Ap,. For each further 
increase in flow rate the increase in 
voidage is such that this constancy of 
Ap, remains intact. This well con- 
trolled condition of equilibrium is 
obviously the requirement of an ideal 
system found only seldom in prac- 
tice. 


Depending on the system and its 
characteristics various types of ir- 
regularities may be experienced. The 
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Figure 1—Pressure drop across a bed 
of ideally fluidizing solids. 
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Figure 2—Pressure drop across a 
slugging bed of fluidizing solids. 





Log G 
Figure 3—Pressure drop across a 
mildly channeling bed of fluidizing 
solids. 


Log ap 


Log G 
Figure 4—Pressure drop across a 
severely channeling bed of fluidizing 
solids. 


more important are slugging and 
channeling. Figure 2 pertains to a 
slugging system. Again normal fluid- 
ization sets in at “a” as the pressure 
drop attains a value of Ap,, but 
constancy prevails only for a brief 
flow interval, and at point “s” the 
pressure drop increases beyond Ap,. 
The increase is caused by wall fric- 
tion of slugs of solids, pushed upward 
by the gas. Since the size (principally 
diameter) of the slugs is of the order 
of magnitude of the vessel diameter 
it is apparent that the phenomenon of 
slugging is found essentially in small 
diameter equipment. 

Slugging is objectionable in process 
studies for several reasons. Probably 
the most important is the effects 
which slugging has on the effective 
space velocity and thereby contact 
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time; the other objections pertain to 
faulty indications as to heat transfer 
and mechanical stability of the solid 
particles. For these reasons artificial 
means such as baffling and other 
modifications in apparatus construc- 
tion have been devised, to eliminate 
slugging. 

The behavior of mildly channeling 
beds is indicated in Figure 3. Bed 
expansion occurs again at point ee 
which, in this instance, corresponds 
to a pressured rop Ap,, in excess of 
A p,- The excess is due to the effort 
that is required to break up the bed. 
Once this is accomplished there ‘is a 
relatively sudden decrease in pressure 
drop (may then go below A p,), 
indicating that a substantial fraction 
of the gas traverses the bed by way 
of regularly formed channels. As 
flows continue to increase the inter- 
nal agitation may become sufficient to 
destroy the channels and thus lead to 
a normal fluidizing bed. If the par- 
ticle characteristics are such as to 
encourage very strong channeling the 
increased fluid velocity may not be 
sufficient to render the bed homo- 
geneous, but may merely be enough 
to rearrange the pattern ‘of the chan- 
nels, leading to sudden breaks in the 
graph, representing the pressure drop- 
flow relationship. This is indicated in 
Figure 4. Such a solid will obviously 
not yield to fluidization under normal 
conditions. 


From a process development point 
of view channeling is very objection- 
able (and possibly more so than 
slugging ), as both space velocities and 
general contacting as well as heat 
transfer charactersitics become great- 
ly obscured. As with slugging we 
may have to rely on artificial means 
to lessen the effect of the phenome- 
non. Unlike slugging channeling may 
occur in both small and large scale 
units. Particle characteristics are a 
major factor in defining and con- 
trolling channeling. From the brief 
consideration of the channeling phe- 
nomenon it follows that incipient 
fluidization in channeling solids re- 
quires higher flow rates than in beds 
of non- channeling solids. 


Having examined the various phe- 
nomena that may possibly occur in 
solid-gas fluidized systems it will be 
instructive to examine briefly their 
relative sequence in the framework 
of the entire fluidized state. Chan- 
neling and slugging are principally 
phenomena that occur in solid-gas 
systems. Solid-liquid systems may 
however channel too, but only under 
the most extraordinary circumstances. 
Particle and general system character- 
istics which are responsible for chan- 
neling in solid-liquid systems (and 


for that matter also in solid-gas sys- 
tems) are only vaguely recognized. 
Solid-liquid systems may also have 
tendencies to slug, but only then 
when the largest density differences 
between solids and liquids are in- 
volved. Thus lead shot fluidized 
water may under certain conditions 
lead to slugging. 

With solid-gas systems if at all 
channeling is most likely to set in 
near or slightly ahead of the point of 
incipient fluidization. Depending on 
particle characteristics mainly, chan- 
neling may endure over a wide range 
of flow rates. Thus mildly channeling 
solids may return to normal fluidiza- 
tion relatively soon, whereas with 
other solids aggregative fluidization 
or even the dispersed phase may re- 
sult after the fixed bed. Obv iously i in 
the latter case the solids do not 
fluidize at all and the pressure drop- 
Huid flow behavior of such a system 
is probably in accord with that given 
in Figure 4. 


Onset of Fluidization 


One of the major problems asso- 
ciated with the design of a fluidized 
system is to predict as closely as pos- 
sible the required mass velocity at 
which fluidization will set in. Since 
in catalytic reaction systems the fluids 
passing through the catalyst undergo 
physical property changes it is im- 
portant to understand how fluid and 
solids physical properties will in- 
fluence the state of fluidization. Re- 
lationships bearing on this important 
question have become available. 


From a consideration of laminar 
flow through granular fixed beds and 
the pressure drop behavior at the 
onset of fluidization there results the 
following expression, permitting es- 
timation of the required mass velocity 
to initially fluidize a charge. Thus: 


0.005 D,2ge pr (ps — pr) be? Ems 
were (1 — Ems) F 
(2) 
In this equation all quantities ex- 
cept shape factor, ¢, and minimum 
fluid voidage, ¢,,¢ are relatively easily 
available from general process specifi- 
cations. In view of this it is of value 
to have observed that for a great 
variety of particles and materials the 
dimensionless product 


9 Ems* 
% 1 — €ms 
may simply be expressed as a function 
of the modified Reynolds number “). 
Thus if this is represented by C, 
Equation (2) takes the form 


C D,? gc pr (ps — pr) 
Ub 


Gini ss 


(3) 
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Evaluation of C, inclusion of g, in 
th experimental constant and ex- 
pressing Dp in inches, pp and ps in 
ounds per cubic foot and wu in centi- 
poises has led to 
> om . \]0.94 
Goi = 688 Dy! [pr ee (4) 
It is to be noted that the dimen- 
sional character of Equation (4) is the 
result of the units chosen and the 
inclusion of g, in the experimental 
constant. The equation is applicable 
to non-vesicular as well as vesicular 
particles and general (tested) limits 
of applicability are as follows: 


D, = 0.205” — 0.002” 

De = 0.5" — 6.0” 

pr = 0.005 — 62.3 Ib./ft.3 
ps = 70-673 lb./ft.3 

w= 0.022 — 2.42 Ib./hr.ft. 


The equation is directly applicable to 
systems where the Rey nolds number 
does not exceed a value of 5.0. For 


DpG 


u 


> 5.0 a correction factor has 


been proposed which extends the 
range of the equation and renders it 
applicable to solid-liquid systems. For 
mixed size solids particle diameter was 
originally defined by 


D, = Gy 4-2 date Te 646 dn Xn (4a) 


where d,, d, ---- d, represent geo- 
metric mean diameters as calculated 
from adjacent sieve openings and x,, 
X) +++. +X, are weight fractions of 
individual size components. Equation 
(4a) is satisfactory for giving ap- 
proximate average diameters of mixed 
size beds, provided the range of sizes 
of the individual components in the 
bed is essentially not in excess of 100 
percent. For more complex mixtures 
it is recommended that particle 
diameter be evaluated according to 


Dh « (4b) 


yr Xn 


da 

The satisfactory agreement be- 
tween values predicted by Equations 
(2) and (4) as pertains to a silica sand 
and a heavy process catalyst is indi- 
~~ in Figure 5. For convenience 
Equation (4) (as applied to gas 
fluidization) has been condensed into 
the form of monograph, Figure 6. 


Example:—A column of 4-in. di- 
ameter is charged with a granular 
catalyst of av erage particle size Dp = 
0.00430 inch to a height of 1.06 ft. It 
is fluidized with air of 80°F. and 734 
mm Hg barometric pressure. The 
catalyst material density ps = 300 
lb./ft.. For these conditions the mini- 
mum fluidization mass velocity can 
be found. 
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Figure 5—Comparison of Equations 


(2) and (4). 


Either Equation (4) or the nomo- 
graph can be used. Fluid density 
28.9 X 734 X 492 


sewed ee = 0.071, and 
359 X 760 x 540 


r= 
(ps — pr) pr = 
(300 — 0.071) X 0.071 = 21.3. 


Gas viscosity u may be taken as 0.018 
centipoises. Using the nomograph this 
vields Gj, = 20 Ib./ft.*hr. 


Bed Expansion During Fluidization 


Whereas the fluid velocity for the 
onset of fluidization may be predicted 
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with considerable accuracy the state 
of bed expansion is not as reliably 
predicted. The reason for this be- 
comes apparent when we recall that 
at the point of initial fluidization 
many systems (except those showing 
severe channeling tendencies) behave 
nearly ideally, whereas at elevated 
flow rates inherent system imperfec- 
tions and deviations from ideal be- 
havior become increasingly _ pro- 
nounced. In line with this we find 
that systems that fluidize particulately 
(solid- -liquid systems, primarily) will 
obey the laws of the fixed bed much 
more readily and over a wider range 
of flow than will the more irregu- 
larly behaving solid-gas systems. 

In order to examine the extent to 
which bed expansions will follow the 
law that applies to fixed beds (Equa- 
tion (2)) it is merely required to plot 
fluid mass velocities versus values of 


e 


(l1—e)’ 
taneous voidage values of the beds, 
associated with any particular value 
of G,, the fluidization mass velocity. 
If the fixed bed law is obeyed then, 
in accord with Equation (2) , such a 
plot should yield a straight line of 


where ¢ represent instan- 


slope—1.0, when logarithmic coor- 
dinates are used. Extensive data 
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Figure 6—Nomograph, Equation (4). 
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analyses have been made along this 
line and the results are indicated in 
Figure 7. Branch A pertains to solid- 
gas and branch B to solid-liquid sys- 
tems. It appears that the gas- fluidized 
systems deviate progressively more 
as particle diameter decreases, where- 
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Figure 7—Deviations from fixed bed 
behavior. 


as the liquid-fluidized systems are in 
substantial agreement with the fixed 
bed flow equation. 

The trend of the data of branch A 
has been the subject of some discus- 
sion), Thus Morse explains the 
deviation from Equation (2) by as- 
suming solids flocculation, that is 
formation of complexes composed of 
individual particles. It is believed that 
it is the complexes which are hydro- 
dynamically active in the charge, 
rather than the individual particles. 


Particle 


The interpretation demands further 
that the size of the complexes in- 
creases with fluid rate. This is some- 
what difficult to visualize, except 
perhaps that in a sense this may be 
looked upon as an increase in the 
state of aggregation. It is highly sug- 
gestive that the difference in expan- 
sion characteristics between gas- and 
l:quid-fluidized systems, as indicated 
by Figure 8 may be associated with 
the characteristics of particulate and 
aggregative bed behavior. 

ndividual expansions of gas-fluid- 
ized systems are readily obtained from 
the data of branch A by way of 
graphical method. The problem is 
simplified by the generalizations of 
by Figure may be associated with 
structed on the basis of a wide range 
of data. For the purpose of estimating 
expansion ratio we calculate first the 
value of G,,- from either Equation 
(4) or the nomograph, Figure 6. 
From Figure 8 we will then obtain 
the value of , the fluidization effi- 
ciency and from Figure 9 we obtain 
the value of the expansion ratio. For 
the purpose of defining the range of 
reliable evaluation of the expansion 
ratio the tested limits and applicability 
of the flow equation are indicated in 
Figures 8 and 9. Fluidization effi- 
ciency has been defined by 

Gr - G 
a 

wherein Gy is the mass tow required 
to fluidize a bed (that is expand the 
bed and keep the particles in motion) 
and G, is a ficticious mass velocity, 
in accord with the flow equation, to 
expand the bed to the same average 


Diometer, inch 
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Figure 8—Fluidization efficiency 


74. 


versus reduced mass velocity. 


voidage as pertains to Gy, but without 
internal particle movement. 

For predicting expansion of liquid- 
fluidized beds Equation (2) will apply 
up to a voidage of about 80-85% 
This is indicated by the data of 
Figure 10 which pertains to small 
particles of sand and catalysts. An 
important ang for applying 
Equation (2) to liquid-fluidized sys- 
tems is _ Ae flow must be laminar, 


Dp G 


that is -< 10. For larger values 


u 

of the Reynolds number (i.e. the 
transition and turbulent range) no 
entirely satisfactory analytical 
methods have been devel oped. Since 
the writing of this paper an impor- 
tant contribution by J. F. Richardson 
and W. N. Zaki, Trans. Inst. Chem. 
Eng’s, London, 32, 35, (1954) has 
come to the attention of the writer. 
Their contribution is most significant 
and useful for the prediction of the 
expansion of liquid-fluidized columns 
up to Reynolds numbers in excess of 
500. 


Example:—For the preceding ex- 
ample the expanded bed height can 
be estimated if the charge is fluidized 
with 100 Ib./ft.2hr. of air. 

From the preceding 
Gur = 20 Ib./ft.2hr.; hence G,/Ginp = 
100/20 = 5. From Figure 8, ny = 0.59 
and from Figure 9, R = 1.20. Hence 
expanded bed height will be 1.06 x 
120) = 4.27 Geet: 


example 


Heat Transfer 


In fluidized systems (just as in 
fixed beds) one must consider the 
problem of heat exchange between 
solids and fluids as well as between 
the bed as a whole and its environ- 
ment. The relation between the two 
phenomena and their respective rates 
are responsible for establishing ther- 
mal equilibrium conditions and tem- 
perature gradients in solid-fluid sys- 
tems. Thus if we consider for 
example an exothermic reaction with 
a considerable heat release on the 
catalyst surface the temperature pat- 
tern through the bed, starting from 
the catalyst surface through the gas 
films to the cooling wall will depend 
on the relative heat removal rates that 
prevail. 

The problem of solids-fluid heat 
transfer in fluidized systems (as well 
as the problem of solids-fluid mass 
transfer) has not yet received suffi- 
cient attention to justify elaboration 
and presentation of proposed results 
with generalization in mind. This is 
probably due to the considerable ex- 
perimental difficulties that are asso- 
ciated with this problem, especially 
the estimation of solids skin. tempera- 
tures. As for heat transfer between 
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bed and environment the difficulties 
are fewer and consequently a fair 
number of studies have been made. 
Observations of Baerg and coworkers 
(6) have indicated that a sharp in- 
crease in heat transfer rates results 
when a fixed bed transforms into the 
fluidized state. As mass velocities in- 
crease coefficients continue to rise, 
though at a diminishing rate and 
finally a maximum coefficient is 
reached. The procession of this phe- 
nomena is indicated in Figure 11. The 
indicated curve picture is typical and 
its analysis is relatively easy, assum- 
ing a certain heat transfer mechanism 
(7), The value of the maximum co- 
efficient will depend primarily on the 
size of the particles, the fluidizing gas 
and its thermal conductivity. Order- 
of-magnitude-wise the maxima fall 
within the range of 50 to about 125 
BTU /hr.°F.ft.?. 

Heat transfer coefficients pertain- 
ing to heat flow (heating and cool- 
ing) between bed and exterior sur- 
faces have been correlated ‘7) and the 
following equation is suggested. 


h D, Cs . D's 0.5 04 
' sae{2h 
k k 
GD" 
( oe ") (5) 


The equation is based on the prin- 
cipal literature data in this field (8, 9, 
10, 11) and experimental limits as per 
data used for arriving at the correla- 
tion are as follows: 


dD, = 0.0015” — 0.0335” 
D, = 2,0” — 4.73” 
gases = air, COs, He, He, A, CHy, He + 
No, Town Gas 
The equation is dimensionally 


homogenous and any consistent set of 


units may be used. 

Limits of flow range for the in- 
dividual data varied from about 20 to 
1500 lb./ft.2hr. Impressive as this may 
appear the limitation as to flow range 
as far as Equation (5) is concerned 
may not be defined wholly satisfac- 
torily by the above mass velocity 
range. This is so because fluidization 
mass velocity must be judged by way 
of particle size. Thus it is more sig- 
nificant to specify a range of reduced 
mass velocities for the equation and 
in line with this it is recommended 
that the equation be limited to appli- 
cations for which G;/G,,, < 20. Ad- 
herence to these high values of re- 
duced mass velocity will lead to 
evaluation of y and R in a range for 
which the curves in Figures 8 and 9 
are dotted (that is they have been 
extrapolated by calculation, and not 
truly established by tests). Neverthe- 
less the calculated heat transfer 
values are closely representative of 
experimental data. 


Approximate Limit of 
Bed Expansion (D,=2.5") 


Particle Diameter, inch 


according to Flow Equation 
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Figure 9—Fluidization efficiency versus expansion ratio. 


Example:—For the conditions of 
the preceding example the possible 
bed-exterior wall coefficient can be 
estimated if c, = 0.2 BITU/°F., lb. 
and the thermal conductivity of the 


air is 0.02 BTU/hr.°F ft. 


(= I dy "\" 
LR 


100 x 0.00430 x 0.59 \°** - 
12 < 0.018 & 2.42 * 1.20 = 0,72 
Cs Ps D,!'5 go's 0.4 s 
k = 
0.2 X 300 X 0.068 x 2.04\°" 
0.02 = 11.1 


and 


hn = O16 X Mt x 0.72 x 12 x 0.02 
a 0.0043 - 


71.5 BTU/hr. °F.ft.2. 

The curve picture of Figure 11 
pertains both to internally as well as 
externally heated or cooled beds, 
However for most beds with internal 
heat transfer elements coefficients 
seem to be considerably higher than 
for beds with external heat exchange 
surfaces. Of course in order to esta- 
blish a definite and valid relationship 
between the two cases the compari- 
son must be made on the same bed, 
that is the unit must have been con- 
structed such that it carried internal 
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and external heat transfer elements 
simultaneously. Moreover tempera- 
ture measurements must have been 
made for both cases. Except for the 
data of Toomey and Johnstone (°) 
such information is not available in 
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the literature. Thus their data have 
been used to relate the ratio of the 
internal to the external heat transfer 
coefficients with the fluidization mass 
velocity (Figure 12). It is interesting 
to observe that as flows increase the 
ratio tends to approach unity. The 
characteristics of the plot point of 
course to significant observations re- 
garding internal gas-solid distribution 
in the bed. More important however 
is the fact that such data may possibly 
provide us with a means of estimating 
the coefficient in (small diameter) 
beds with internal surfaces once the 
coefficients for beds with external 
elements can be estimated indepen- 
dently. 

Probaly the most important study 
of heat exchange between an im- 
mersed heat transfer element and a 
surruunding fluidized bed was made 
by Baerg et al ‘®). Order of magni- 
tude-wise the coefficients are higher 
than coefficients pertaining to bed- 
exterior wall systems, operating under 
comparable flow rates. This is clearly 
in line with the data of Figure 12. 
Another study of fluidized bed heat 
transfer, carrying an internal ex- 
changer surface was made by Vree- 
denberg ‘'?), His data are of particu- 
lar interest since he could demonstrate 
that the position of the immersed 
element in the bed has a profound 
bearing on the value of the heat 
transfer coefficient. In other words 
the Vreedenberg study emphasizes 
clearly that an improved understand- 
ing of fluidized bed heat transfer can 
only result from a better under- 
standing of the internal mechanism 
of the fluidized state. 


Ratio of Internal to External Coefficient 
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Figure 12—Internal-external fluidized 


bed heat 


transfer coefficient 
mass velocity. 


versus 


Nomenclature 


Cs 
D, 
D, 
& 
G, 
G; 


Gm 


heat capacity of solids. 

function of Reynolds number. 

particle diameter (inches in Equa- 

tion 4). 

vessel diameter. 

acceleration due to gravity. 

ficticious mass velocity, required to 

expand but not fluidize a bed. 

mass velocity required to fluidize a 

bed. 

mass velocity required to incip- 

iently fluidize a bed. (Ib./ft.*hr. in 

Equation 4). 

coefficient of heat 

wall. 

gas thermal conductivity. 

bed height. 

slope, when log G; is plotted versus 
(i =) 

log 


transfer, bed- 


é&3 


= bed expansion ratio. 


pressure drop. 
Subscript t—theoretical, c 
neling, s—slugging. 


chan- 


€ = bed voidage. 

€mi = minimum fluid bed voidage. 

7 = fluidization efficiency. 

KM == fluid viscosity. 

pr = fluid density. 

_(Ib./ft.° in Equation 4), 
ps = solids density. 

@: = particle shape factor. 
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Wetted Wall Packing — 


A New Distillation Column Packing’ 


A. R. BANCROFT? and H. K. RAE* 


The performance of a wetted-wall type of packing consisting of 
closely spaced parallel vertical sheets of aluminum has been investigated. 
The liquid flows down the surfaces of the sheets as a thin film in contact 


with the rising vapor. 


Tests were conducted in a 6 in. x 6 in. square 


column 50 ft. high at total reflux with the systems H,.O-HDO and ethanol- 


isopropanol. 
used. 


Four sheet spacings ranging from 3/16 to 1/2 in. were 


Throughputs as high as 2600 lb./hr. sq. ft. (vapor velocity of 20 
ft./sec.) can be obtained with a pressure drop of less than 0.10 in. of 


water per ft. of column height. 


HTU values range from 2 to 8 ft. The 


pressure drop per transfer unit is lower than for any other type of 
distillation column at the same throughput, although many other pack- 


ings give much lower HTU values. 


The maximum flow rate for this 


packing is somewhat greater than for plate distillation columns and is 
much greater than for most other packings. 


The data are also compared with correlations for other types of 


wetted-wall columns. 


Introduction 
oe paper describes the operating 
characteristics of a new type of 
distillation column packing working 
on the wetted-wall principle. The two 
important characteristics of such a 
packing are a very low pressure drop 
and a high flow capacity. The devel- 
opment of this packing was under- 
taken in connection with studies of 
the economic production of heavy 
water by the distillation of natural 
water. 

Benedict ‘!) in his excellent review 
of heavy water production technology 
discussed water distillation and point- 
ed out some of the desirable column 
characteristics that would result in 
low energy and capital costs for the 
process. Selak and Finke?) have re- 
ported the operation of the Manhat- 
tan Project distillation plants and in- 
dicated how the high costs which 
were encountered might be reduced. 
In an earlier paper Rae‘) discussed 
the production of heavy water and 
has shown under what conditions 
water distillation could be attractive. 

Both the low value of the relative 
volatility and the low concentration 


of deuterium in natural water result 
in extremely large boil-up rates in the 
primary columns per pound of pro- 
duct and also in a large number of 
theoretical plates. Low energy costs 
can be obtained by vapor compression 
to give efficient heat recovery ‘*) or 
by using low pressure by- -product 
steam from a power plant either 
directly in the columns‘) or in the 
reboilers"!). In each case it is desir- 
able to have the lowest possible pres- 
sure drop in the columns, and a figure 
of 0.1 in. of water per theoretical plate 
is suggested for the vapor compres- 
sion case‘), High flow capacity and 
low heights equivalent to a theoretical 
plate are required to give a reason- 
ably small total column volume. 


An examination of the pressure 
drop data for various distillation 
columns (3) indicated that the wetted- 
wall column was the only type with 
the desired operating characteristics. 
Experimental data have been obtained 
by a number of inv estigators for 
small-diameter single-tube (6.7.8.9) and 
for small-diameter multiple-tube 
(10,11) wetted-wall columns. Scale-up 
of the wetted-wall principle by using 


1 Manuscript received May 13, 1957. 
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closely spaced parallel vertical sheets 
was envisioned by the authors to be 
more practical than parallel or con- 
centric tubes. 

Estimates for heavy water produc- 
tion using the wetted-wall column in 
a vapor compression cycle and basing 
the design of the parallel vertical sheet 
packing on cerrelations of single-tube 
data, gave a cost of about $30 per 
pound (!?), At that time this figure 
was sufficiently attractive to initiate 
the development of a suitable wetted- 
wall packing, but more recently the 
USAEC have announced a price of 
$28 per pound*. However, the devel- 
opment work was continued because 
of the possibility of alternative appli- 
cations for the packing. 


Single-tube Wetted-wall Columns 

In the streamline flow region for 
smooth wetted-wall channels West- 
haver ‘!3) has shown that the HTU 
should be given by 


(HTU )y = 


0.1215 Re)y (Sc)y + — (1 
-1Z19 Th Py CSC by 
fad SN SE te ek : 


The second term in the bracket is 
negligible if (Re), is greater than 100. 
In the turbulent flow range it has 
been found (67) -that the Chilton- 
Colburn analogy for mass, heat and 
momentum transfer (14) will correlate 
the data for small-diameter single-tube 
wetted-wall distillation columns. Thus, 
f, defined by the usual Fanning fric- 
tion equation (!5), and jp are related 
by 
= 0.033 (Re),~°:?4 (2) 
Rearranging, the HTU is given by 
(HTU) 
All the available data indicate that for 
distillation in these columns (HTU),, 
is equal to (HTU), and the process 
is controlled by diffusion in the vapor 
phase. In the few cases where it has 
been possible to obtain values of 
(HTU), ‘“8) these have been very 


©The process used by the USAEC is based on 
the chemical exchange of deuterium between 
water and hydrogen sulphide (1). 


= 30.3 rn (Re)v®28 (Sc)y%8? (3) 
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Figure 1—Diagram of distillation column. 
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small, indicating negligible resistance 
to mass transfer in the liquid phase in 
these wetted-wall distillation columns. 

If the definitions of f and jp are 
substituted in equation (2) and the 
pressure drop is expressed in inches 
of water, then by rearranging, the 
pressure drop per transfer unit is 
found to be 


G2 (Sc)? 3 
APry = (4) 


5.2 gp, 

From equation (4) it will be seen 
that very low values of the pressure 
drop per transfer unit can be achieved 
at high vapor rates. For the HXO-HDO 
system a pressure drop of 0.12 in. of 
water per transfer unit is obtained 
at a vapor rate of 3600 lb./hr. sq. ft. 
Applying equation (3), the (HTU),, 
under these conditions would be 1.6 
ft. for a column diameter of 0.5 in. 
Such values of AP;y, G, and (HTU),, 
are typical of the desired operating 
characteristics of a large scale wetted- 
wall column for water distillation. 


Early Development 

Three requirements for the efficient 
operation of a large-scale wetted-wall 
column are: 

(i) obtain as high a ratio of wetted 
perimeter to cross-sectional area as 
possible; 

(ii) distribute the liquid uniformly 
over all of this wetted perimeter; 

(iii) maintain this uniform distri- 
bution of the liquid film throughout 
the length of the column. 

The use of closely spaced parallel 
vertical sheets was considered the 
most economical arrangement to yield 
the desired small hydraulic radius. 

Initial development work was done 
with small-diameter columns where it 
was found that flat imperforated 
sheets of aluminum gave much better 
performance than fine-mesh screen 
both as to uniformity of wetting and 
low pressure drop. Various configura- 
tions of parallel vertical sheets were 
tested in 6” x 6” and 12” x 12” square 
lucite columns using air and water to 
study liquid distribution, pressure 
drop and flooding phenomena (‘!®), A 
simple and structurally strong pack- 
ing unit was developed and used in 
the distillation work described in this 
paper. 

Operation of the distillation column 
was planned to investigate the vari- 
ation of the pressure drop and. the 
height of a transfer unit with flow 
rate and sheet spacing, and the effect 
of length of the individual packing 
units on column performance. Initi- 
ally, only water distillation runs were 
planned, but the column was modified 
to operate with an alcohol system as 
well when the column performance 
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Figure 2—Column head. 


did not agree with existing correla- 
tions and the effect of varying the 
physical properties of the liquid and 
vapor became of interest. 


Description of Apparatus 

The distillation column is 6 in. sq. 
and has an effective packed height of 
47.3 ft. Vapor is generated in a steam- 
heated reboiler, passes through 2 
packing, and is totally condensed 
the column head. The condensate i 
preheated before being returned to 
the column. The system operates at 
total reflux; no provisions are made 
for adding feed or withdrawing pro- 
duct. A schematic diagram of the 
apparatus is given in Figure 1 . Photo- 
graphs of the column base and top are 
shown in Figures 2 and 3. 


Packing 

The packing units have dimensions 
6” x 6” x 24”. A typical unit is shown 
in the photograph, Figure 4. Each unit 
is assembled from a number of alu- 
minum sheets (1/16” x 6” x 24”) and 
six “spacers” (1/16” x 1” x 6”). Figure 
5 shows the manner in which the 
spacers are fitted into lugs on the 
sheets to form a rigid unit. The spac- 
ing, which ranges from 3/16 to 3 in., 
is measured between adjacent faces 
of the sheets, and the number of 
sheets per packing unit are as follows: 
Spacing—inches 3/16 4+ #§ 34 
Sheets per packing unit23_ 19 13 11 


Column 

The column shell is composed of a 
number of 4-ft. sections of 4-in. alu- 
minum, 6-in. sq. in cross section, 
mounted in guides which allow for 
column expansion and permit ver- 
tical alignment. Two-in. magnesia 
block insulation is used to cover the 
entire column. The calculated heat 
loss is less than 1% of the total heat 
transport up the column at the lowest 





flow rate. The column sections are 
joined and sealed by means of a bell- 
and-spigot-type joint filled with Aral- 
dite 101 resin. This resin makes a 
satisfactory seal provided the column 
pressure drop does not exceed 12 in. 
of water, and can be easily removed 
from the joint for dismantling the 
column to change the packing. At the 
column base and top, where vapor 
enters and leaves the packing, the 
shell cross section is enlarged to per- 
mit the vapor to flow through the 
sides of the packing units instead 

the end. This method causes a mini- 
mum disturbance to the liquid flows. 


Figure 3—Column base. 


Distributor 

Liquid is distributed evenly to the 
top of a packing unit extending above 
the vapor outlet section of the ‘column 
by means of a box-type distributor. 
In the bottom face of the box holes 
(0.030-0.040"D) are drilled on a square 
grid to match the spacing of the 
sheets, giving a line of holes spaced 
closely together (} to 3”) above each 
sheet. Larger holes on the outer row 
improve the unifor mity of distribution. 
The distributor is raised slightly above 
the packing unit (about 1/64”) to 
prevent the packing sheets from block- 
ing the holes. 


Operation 

The flow r rate in the column is re- 
gulated by a Taylor Fulscope Con- 
troller which is sensitive to pressure 
at the bottom of the column (and, 
therefore, flow) and controls the 
flow of steam to the reboiler. Column 
pressure, generally in the range 0-10 
in. of water, is transmitted from the 
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base through liquid seals to the in- 
strument. The condensate flows by 
gravity to a covered 20-gallon alu- 
minum tank, from which it is pump- 
ed through a rotameter and preheater 
to the distributor at the top of the 
column. Since the system operates at 
total reflux, the reflux rate is kept 
constant by means of a float-operated 
valve which maintains a constant level 
in the tank. The reflux temperature 
is regulated by a self-acting controller 
in the liquid line, which controls the 
steam flow to the preheater. 

Stainless steel or aluminum is used 
for all the equipment to minimize any 
corrosion and provide an absolutely 
clean system. 

Water cooling at various points 
minimizes the loss of vapor from 
vented portions of the system. 

The column flow rate is measured 
as liquid reflux by means of the rota- 
meter, corrected for a small amount 
of liquid which is returned to the 
tank from the vapor disengaging sec- 
tion. This is reported as the mass flow 
rate of the vapor and can be taken as 
constant throughout the column and 
equal to the mass flow rate of the 
liquid for the two systems studied. 

The column pressure is measured 
at five different elevations. Pressure 
taps are 3-in. tubing flush welded to 
the column shell at a 45° slope, and 
are connected to simple Meriam mano- 
meters filled with S.G.= 1.0 liquid. 
Pressure readings are taken to the 
nearest 0.05 in. Sy stem temperatures 





Figure 4—Packing unit, 4 in. spac- 
ing. 
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Figure 5—Assembly of packing unit. 


are recorded for operating conveni- 
ence. Liquid samples are taken for 
analysis from the column top and the 
reboiler. 


All mass transfer and pressure drop 
measurements are taken during 6-hr. 
tests. Within an hour of start-up the 
system reaches steady flow rate and 
liquid level conditions, and within 3 
hr. the top and bottom concentrations 
reach steady state. 


Results 
Liquid Distribution 

Although it was not possible to 
measure the distribution of the liquid 
on the packing units in the distilla- 
tion column, the results obtained 
from measurements at room temper- 
ature in air are of interest because of 
the importance of uniform distribu- 
tion in the efficient operation of the 
column. 


For the water system, aluminum 
etched in a solution of sodium hy- 
droxide (approximately 5%) present- 
ed an easily-wetted surface which 
minimized any channelling of the 
liquid and promoted complete wet- 
ting of the surface. 


The action of the distributor was 
observed to result in the formation 
of separate streams under each hole, 
but these streams united to form a 
continuous film covering the width 
of the sheet within an inch or so 
from the top. This film was found 
not to rupture or to channel for a 
distance of several feet at the flow 
rates used in the distillation column. 
To reduce the possibilities of chan- 
nelling on very long sheets a method 
of redistributing the liquid at inter- 
vals was used which consisted of 
placing the sheets of each packing 
unit at right angles to the sheets of 
the unit above it. (Refer to Figure 4). 
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Using ten units stacked in this man- 
ner with a distributor at the top 
satisfactory redistribution at the bot- 
tom of each unit was observed in a 
number of out-of-column tests. Feed- 
ing water to a single sheet of the top 
packing unit resulted in a gradual 
spreading of the water to cover as 
many as seven sheets of the bottom 
unit. This indicated that the redistri- 
bution results in some mixing of the 
liquid from different sheets and should 
tend to correct any non-uniformity 
in the flow pattern. 

The liquid flow rate from the bot- 
tom of each sheet of a packing unit 
was measured using a suitable collec- 
tion trough, the packing unit being 
placed below a second unit to which 
the water was fed from a box-type 
distributor. A typical distribution 
diagram is shown in Figure 6 for a 
sheet spacing a 3/16 of an in. The 
distribution for other spacings was 
similar. The liquid flow rates increas- 
ed slightly from one side of the unit 
to the other, and the average devi- 
ation from the mean flow was less 
than + 10%. Very little water flowed 
on the two outside sheets at one side 
and on one sheet at the other for this 
particular unit. However, at this stage 
of the work the distributor was not 
properly fitted to the packing, and 
since it was felt that this situation 
could be easily improved no addi- 
tional tests were made. 

In the distillation column the tur- 
bulence of the rising vapor and the 
higher temperature of operation were 
expected to improve the uniformity 
of distribution of the liquid over the 
results reported here. 


Entrainment and Flooding 

When liquid flowed as a film down 
a plane vertical sheet it accumulated 
at the bottom of the sheet until there 
was enough to form drops. If another 
sheet was placed below the first and 
touching it at right angles, then the 
liquid drained onto the second sheet 
rather than dropping off. When there 
were closely spaced “drains” at the 
bottom edge of the sheet the accumu- 
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lation of liquid was small, but when 
the “drains” were widely spaced the 
accumulation was greater. This de- 
scribes the liquid behaviour at the 
contact plane of two packing units. 

For 3-in. spaced packing consider- 
able liquid accumulated at the higher 
flow rates and was easily blown off 
the sheets by the rising vapor. The 
small droplets thus formed were car- 
ried upward with the vapor until a 
collided with and were absorbed i 
the liquid film. At lower sheet eis 
ings where less liquid accumulated 
droplets were not so readily blown 
from the bottom of the sheets, and 
entrainment did not begin until higher 
liquid and vapor rates had been 
reached. Flooding occurred when the 
liquid bridged between the sheets 
at the contact plane of two packing 
units and the vapor bubbled through 
the liquid. The transition from en- 
trainment to flooding was gradual in 
most cases. 

With the aqueous system droplets 
were entrained at lower flows than 
flooding for all spacings, and no dis- 
continuities were found in the column 
performance as the flow rate was in- 
creased. However, the alcohol system 
behaved differently at both 3- and 
3/16-in. spacing. At the 4-in. spacing 
as the flow rate was increased small 
droplets were entrained and carried 
about 8 in. up the sheets. When the 
flow rate reached a definite value 
(2100 Ib./hr. sq. ft.) only large drop- 
lets were formed and, because of 
their increased size, were carried up 
only an inch or so. This change in 
entrainment was reflected by “breaks” 
in both the pressure drop and HTU 
curves (Figures 8 and 9). For the 
3/16-in. spacing flooding occurred at 
a definite flow rate (1750 lb./hr. sq. 
ft.) and no entrainment was observed 
at lower flows. This flooding condi- 


tion started at the bottom packing 
units and a 10% increase in flow was 
required to flood the entire 48 ft. of 
packing. This is probably due to tur- 
bulence created at the vapor entrance 
to the column permitting flooding to 





12 13 14 15 16 17 18 1920 21 22 23 


SHEET NUMBER 


Figure 6—Distribution of liquid flow on sheets of packing unit. Total flow— 
2200 Ib./hr. sq. ft. 
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occur more easily here. As the flow 
rare was increased past this stage the 
pressure drop increased rapidly in a 
continuous manner. 

The maximum flow rate for a pack- 
ing is generally determined by flood- 
ing or entrainment which result in 
an excessive pressure drop or an in- 
creased HTU value or both. Thus, 
this limit can be evaluated by examin- 
ing the HTU and pressure drop plots 
against flow rate (Figures 7, 8 and 9). 
However, the values picked in this 
way are influenced by the shapes of 
the curves, and where changes i in slope 
are not abrupt only an approximate 
figure can be obtained. To standardize 
the ‘ ‘flooding limit” a value of 0.4 in. 
of water pressure drop per transfer 
unit has been set. The flow rates 
which satisfy this condition are listed 
in Table 1. 

The values of the vapor velocity 
and Reynolds number are based on 
the actual cross section of the pack- 
ing available to the vapor. However, 
the values were calculated neglecting 
the thickness (0.0002-0.0005 ft.) and 
the velocity (up to 0.5 ft. per sec.) 
of the liquid films. An error in Re, 
as great as 5% may thus be intro- 
duced, but the effects tend to cancel 
each other and the error is likely 
smaller than this. 


For the aqueous system the flood- 
ing vapor velocity was almost inde- 
pendent of spacing. An average value 
of 22 ft./sec. could be taken as the 
limiting velocity. However, the maxi- 
mum throughput based on the empty- 
column cross section varied from 
2180 to 2720 lb./hr. sq. ft. with pack- 
ing spacing, and the vapor Reynolds 
number changed three-fold from 2970 
to 8430. The flow rates for the alco- 
hol system were slightly lower, but 
the Re, values were similar. 


A limiting flow rate calculated from 
a correlation given by Kamei(!7) for 
open-tube wetted-w all columns was 
found to be about 50% higher than 
that obtained in this packing. The 
difference is probably accounted for 
in the constrictions and irregularities 
at the contact plane between every 
pair of packing units, and also by the 
fact that entrainment, as well as flood- 
ing determined the limiting flow 
rates. 

The first tests in the column, on 
|-in. spaced packing were made using 
a column base where the vapor entered 
the bottom of a packing unit and 
special provision was made for col- 
lecting the liquid immediately below 
the packing to prevent re-entrain- 
ment. This was found to cause pre- 
mature flooding and an_ increased 
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Figure 7—Measured pressure drop for water distillation (sheet spacing shown 
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Figure 8—Measured pressure drop for alcohol distillation (sheet spacing shown 





TABLE 1 
LIMITING FLow RATES IN WETTED-WALL PACKING 
































Sheet Flow Rate Vapor 
Spacing Go Re Velocity 
in. lb./hr. sq. ft. ft./sec. 


H,O0 — HDO system 


36 2180 2970 


3 21.4 

4 2350 2960 an 9 

% 2720 6490 23.4 

6 720 8430 22.8 
ethanol-isopropanol system 

3% 1770 3070 5.7 

4 2540 8750 7.0 











600 1000 1400 1800 2200 2600 3000 
FLOW RATE Ib / hr sq ft 


as third parameter). 




















eee eee ee ie ee ee 





_—+ a 


600 1000 1400 1800 ~ 2200 2600 3000 
FLOW RATE ib / Net sq ft 


as third parameter ). 








18 to pressure drop compared with the 
\ upper sections of the column. As a 
\ result the flooding rate for this spac- 


\ / ing given in Table 1 may be slightly 


low and the pressure-drop curve in 


14 | = La, } ne Figure 7 may be somewhat higher 

\ | | iH than it should be. In addition, the 
| \ | = | HTU curve in Figure 9 may be too 
| oh al | far to the left. 








| | > 26 
L 4 / Pressure Drop 








| 7 The column pressures measured at 
+ A : ; 

4 dh A \ & -- five different locations were plotted 
| oy —e 4— against height from the top and the 

ee ee best straight line was drawn through 
lo | z_. these points. The slope of the line 
T ' we ae eT was reported as the pressure drop per 
a3 WATER | Ho foot of packing. In Figures 7 and 8 


Hs td — Bit. p—— 2 — =o | | the curves for pressure drop against 
veal as Vg 5 | flow rate are shown for both the 

| | water and alcohol systems at different 

chs. ee —————— sheet spacings. In each Figure the 
pressure drop calculated for one spac- 
ing using the Fanning equation for 


$00 1000 1400 1800 2200 2600 3000 turbulent flow in smooth channels (15) 
is shown. The calculated pressure 


FLOW RATE Ib/hr sq ft drop curves for other spacings differ 
Figure 9—Variation of HTU with flow for water and alcohol distillation (sheet from the experimental curves by 
spacing shown as third parameter ). similar amounts. These calculations 
were made neglecting the thickness 
and velocity of the liquid films and 
considered ‘only a straight uniform 
channel; they were corrected for the 
difference in density between air and 
water or alcohol vapor. This correc- 
tion was necessary because the mano- 
meters were zeroed with air in the 
column at room temperature. For 
water distillation the correction re- 
sulted in negative values of the cal- 
culated pressure drop at low flow 
rates, while for the alcohol system 
the correction was positive. 
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At low flow rates the measured 
and “theoretical” values are almost the 
same; at higher flows the measured 
values are consider ably greater. The 
difference can be explained by the 
obstruction to flow presented by the 
edges of the sheets at 2-ft. intervals 
and by the water accumulated at these 
points which reduces the cross-sec- 
tional area of the channels. This resis- 
tance to flow can be calculated for 
dry sheets by considering the number 
and thickness of the sheets, but the 
problem becomes very complicated 
when it is necessary to consider also 
the liquid accumulated on the bottom 
edges of the sheets. The presence of 
the liquid film covering the full length 
of the sheets also causes the measured 
pressure drop to be higher than the 


500 


300 


200 


100 sce calculated value, but this effect is 
1000 2000 4000 6000 8000 relatively small. The large divergence 
(Re)v of the two curves at very high flow 


Figure 10—Comparison of mass transfer data with Chilton-Colburn analogy rates % the result of the onset of 
(sheet spacing shown as third parameter). entrainment and ‘flooding. 
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Mass Transfer 

Samples of tops and bottoms were 
taken every hour from the system 
ope rating at total reflux. For the 
aqueous system the D,O content was 
0.1-0.2 mole %; for the alcohol sys- 
tem the isopropanol content was 25-75 
mole %. Infra-red spectrometry was 
used fox analysis in both systems giv: 
ing an accuracy better than + 0.5% 


The number of transfer units was 
calculated using the Chilton-Colburn 
equation (18); 


1+ (NTU)ow = In = - + 


The relative volatility for HyO-HDO 
separation (at low heavy water con- 
centrations more than 99.9% of the 
deuterium is present as HDO(!®)) 
1.025 (19.20), Based on vapor pressure 
measurements?!) and vapor-liquid 
equilibrium studies ‘??) of the ethanol- 
isopropanol system the relative vola- 
tility was taken as 1.158 at atmospheric 
pressure. The mixture in the column 
normally contained a small amount of 
water which decreases the relative 
volatility of the two alcohols and 
causes a systematic error in the ana- 
lysis. Studies with known amounts 
of water in a small laboratory column 
having about the same number of 
transfer units as the large column 
enabled correction factors to be 
worked out. The results for the 3/16- 
in. spacing were corrected in this way 
but those for the }-in. spacing were 
not because the water content was not 
measured in each run. The resulting 
HTU values at this spacing are about 
ten per cent too high. For both sys- 
tems the relative volatility is constant 
throughout the column. 

The values of the HTU are plotted 
against flow rate for the various spac- 
ings studied with both the water and 
alcohol systems in Figure 9. The 
HTU values tend to increase at both 
ends of the flow range. At low flow 
rates this increase is probably due to 
less uniform liquid distribution, and 
in the case of the alcohol system with 
3/16-in. spacing very poor wetting is 
indicated. At high flow rates the in- 
crease results from the effects of en- 
trainment and flooding. 


The limitations of the column and 
the control system with regard to 
pressure drop prevented low flows 
being reached with the }-in. spacing 
and high flows with the 3/16-in. 
spacing. 

The HTU data are compared with 
the correlations for single-tube wetted- 


TABLE 2 
PHYSICAL PROPERTIES 





Property H,0 D.O Ethanol Isopropanol 

at 100°C at 100° ~ at 80°C at 80°C 

B. Pt., °C at 1 atm. | “100 i 401. 42 : 78.5 . on 82 3 

p., Ib/cu. ft. 0.0373 0.042 0.100 0.131 

pi, Ib/cu. ft. 59.7 66.3 46.1 45.9 

My, Cp. 0.0125 0.0128 0.0104 0.0091 

My, cp. 0.283 s329 0.42 0.53 

H.O — HDO ethanol-isopropanol 

D, ft?/hr. ¥3 0.2 

D, ft?/hr. 4.0xit- 1.2x10- 

Scy 0.63 0.77 

Sci 28 208 

a 1.025 1.158 

Notes: 


Pv, pi, My and mw; for water, ethanol and isopropanol are taken from the handbooks. 


For D.O p, is calculated from HO vapour density and the ratio of molecular weights, 
p: and my are taken from a report by Cooper 2), and py is extrapolated from results of 


Bonilla, Wang, and Weiner‘? 


Dy, is calculated by the inline of Gilliland’?® and Fair and Lerner“®) in the absence 


of experimental values. 


D, for H2O0 —HDO is taken from data of Wang, 


Robinson and Edelman‘); for 


ethanol-isopropanol it is calculated by method of Wilke‘** 


wall distillation columns (equations 
(1) and (3)) in Figure 10. The HTU 
values are all much greater than 
would be expected from these corre- 
lations. Values for the physical pro- 
perties used in the calculations for 
the two systems are given in Table 2. 
A number of explanations for these 
differences are advanced below but 
it is not possible to specify which are 
most important on the basis of the 
data available. Future work is planned 
which may clarify this situation. 

(1) The diffusion coefficients used 

determine the Chilton-Col- 
burn values may not be accur- 
ate. Information on experimen- 
tal values of D, is sparse, and 
equations for calculating it are 
of limited accuracy. 

(2) The possibility of back-mixing 
of the liquid at the contact 
planes of the packing units may 
increase the HTU values some- 
what. 

(3) Although the liquid distribu- 
tion appears to be uniform, 
liquid channelling and unequal 


flow on each sheet may be 
occurring. 
(4) A large resistance to mass 


transfer may exist in the liquid 
phase with this packing al- 
though other wetted-wall dis- 
tillation columns do not show 
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this. There are no theoretical 
or semi-empirical methods for 
predicting mass transfer in the 
liquid film under distillation 


conditions which could be 
compared with actual perform- 
ance. 


The difference between experimen- 
tal and predicted HTU values is much 
greater for the alcohol system than 
for the water system. This could be 
due to a number of factors. The al- 
cohols have a much lower surface 
tension than water which could give 
rise to different liquid flow patterns. 
The entrainment and flooding be- 
haviour is different for the two sys- 
tems. The rate of diffusion in the 
liquid phase is much lower for the 
alcohol system than for the water 
system (see Table 2) so that any 
liquid- phase resistance to mass trans- 
fer will be greater for the former. 


An interesting comparison with 
commercially available packings can 
be made for the pressure drop per 
transfer unit based on equation (4). 
In ee 11 APyy is plotted against 


92(Sc)?/8 

5 ab Spy 
‘alia and for the wetted-wall 
packing. The conditions for each 
curve are given in Table 3. For most 
of the systems included and the two 
studied in the work reported here the 


for eight commercial 
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TABLE 3 


PRESSURE DROP PER TRANSFER UNIT IN PACKED COLUMNS 
Summary of Data Plotted in Figure 11 


System (Sc) Packing 
H.O-air 0.60 carbon grids, 

Vora 

34” pitch 

cre yssed 

serrated 
H.O-air 0.60 wood grids, 

uy" x1" 

1” pitch, 

crossed, plain 
H.O-air 0.60 ceramic rings 

3” x 3” 

3 gm walls 
methylcyclo- 1.06 Cannon 
hexane-n- Protruded 
heptane packing, 

yy" x \Y" 
ethylene 1.18 Stedman 
dichloride 
benzene 
H,O-HDO 0.63 Spraypak 
NH,-air 0.61 Spiral tile, 

3" x 3” 
NH;-air 0.61 Raschig rings 


” x ” 


Equivalent 


Hydraulic Reference Symbol 


Radius 
ft. 

0.0275 29 

0.030 30 

0.050 30 

0.0043 on” 
aa” 
33* 

0.020 34 

0.030 34 


* The data for this case are reported in terms of height of an equivalent theoretical plate, 


so that the pressure drop used in Fgure 11 is that per theoretical plate. 
there is, however, no very great difference between the H.E 


HTU and the height equivalent to a 
theoretical plate (HETP) are ap- 
proximately the same so that Figure 
11 can also be considered as a plot 
of the pressure drop per theoretical 
plate. 

It will be noted that all packings 
except grids have pressure drops very 
much higher than true wetted- wall 
columns and this is because the vapor 
channels are not uniform and bubbl- 
ing through the liquid occurs. For the 
wetted-wall packing both the pressure 
drop and the HTU are greater than 
in simple wetted-wall columns so 
that the AP, is several times greater 
than predicted by Equation (4). In 
the case of the alcohol. runs the dis- 
crepancy is very large because of the 
much higher HTU values. 

For water distillation and spacings 
of }-in. and greater the packing has 
a lower pressure drop per transfer 
unit at a higher flow rate than any 
other type of packing. The only other 
packing which comes close in this 
connection is carbon grids and the 
wetted-wall packing has a_ pressure 
drop which is 50% lower at a given 
flow rate and can operate at higher 
flow rates. 
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In this system 


y.T.P. and the H.T.U. 


Applications 

The recent reduction in price of 
heavy water to $28/lb. has changed 
the outlook on this development pro- 
ject. Estimates have now shown that 
heavy water cannot be produced eco- 
nomically by using water distillation 
in a wetted-wall column, although 
the process may be useful for final 
up-grading of partially enriched 
water. 

The wetted-wall type of packing 
should offer advantages for vacuum 
distillation systems w here a very low 
pressure drop is desirable, particular- 
ly if a large number of transfer units 
are required; or for other mass trans- 
fer operations where a low pressure 
drop is an advantage, especially when 
coupled with a high throughput 
which reduces column size and cost. 
More development work, however, is 
required before the packing can be 
used in large-scale columns. 

An approximate cost for the pack- 
ing units fabricated here is $35 per 
cu. ft. This cost includes both ma- 
terials and labor. The labor cost can 
be reduced by mass producing the 
pieces, and the materials cost can be 
reduced by using thinner packing 


sheets. It is thought the price of 
packing produced in large quantities 
would be less than $20 per cu. ft. for 
t-in. spacing. 

In connection with packing costs, 
the loading strength of 6” x 6” x 24” 
packing units was determined by test- 
ing the units inside a section of 
column shell. The load limit was taken 
as the load required to buckle two 
sheets, and was calculated in feet of 
packing since the bottom unit sup- 
ports all those above. For 5/16-in. 
spacing the following values were de- 
termined. 


Sheet 0.0625 0.036 0.032 0.025 
Thickness in. 
Load Limit 123 64 58 44 


ft. of dry packing 

The results for other spacings were 
very similar. The amount of alu- 
minum required for packing in a 
column could be kept to a minimum 
by using thicker sheets for the bot- 
tom units and thinner sheets near the 
top. 


Conclusions 

The wetted-wall packing described 
in this paper has a very low pressure 
drop at very high flow rates. For ex- 
ample for water distillation, when the 
spacing between the vertical sheets 
of the packing was § of an in. the 
pressure drop per transfer unit at a 
flow rate of 2600 Ib. per hr. sq. ft. was 
0.3 in. of water. This pressure drop, 
especially at such a high flow rate, 
is very much lower than for any 
other type of packing or any plate- 
type of distillation column. Some grid 
packings have pressure drops which 
are not very much higher than the 
wetted-wall packing but they cannot 
operate at such high throughputs. The 
cost of the packing is expected to be 
less than $20 per cu. ft. and it should 
prove useful for mass transfer ope- 
rations where a low pressure drop is 
required. 


Both the pressure drop and the 
height of a transfer unit are higher 
than would be predicted from cor- 
relations for small-diameter single- 
tube wetted-wall columns. The rea- 
sons for these discrepancies are not 
fully understood, but further devel- 
opment work is expected to help in 
the interpretation of the data. 


The performance of the packing 
has been studied for the H»O-HDO 
system and the ethanol-isopropanol 
system. In the latter case the column 
operation was much less efficient, but 
it has not been possible to explain 
why this should be so. 
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‘Lecognizing the importance of uni- 
form liquid distribution throughout 
th packing many  out-of-column 
measurements were made using water. 
These indicated that a_ reasonably 
umform distribution should be ob- 
tained in the column. 


Nemenclature 

A = cross-sectional area of chan- 
nels 

Ag = cross-sectional area of vertical 
sheets 

D = diffusivity 

d = diameter 

f = Fanning friction factor = 
2 gdF 
4u? 

F = frictional pressure drop per 


foot of channel in feet of fluid 
flowing 


g = acceleration due to gravity 
G = vapor flow rate per unit cross- 
sectional area of channel = 
upy 
G. = vapor flow rate based on 
empty cross-section = 
A 
u py 
As 


(HTU). = height of a transfer unit based 
on over-all resistance in vapor 
phase concentration units 

jp = j factor for mass transfer = 
rp (Sc)?/3 

(HTU), 
(NTU )oy = number of transfer units in 


the column based on over-all 
gas phase resistance 


AP = pressure drop (in. of water 
per ft.) 

APru = (AP) (HTU), pressure drop 
per transfer unit in in. of 
water 
4 rh up 

(Re), = , Reynolds number 

bh 

Th = d/4 or s/2, hydraulic radius 

- = spacing between sheets 

: yr « : ‘ 

Sc = , Schmidt Group 
Pp 

u = mean vapor velocity in chan- 
nels 

y = mole %, more volatile com- 
ponent 

a = relative volatility 

be = Viscosity 

p = density 

Subscripts 

1 = top of column 

2 = reboiler 

\V = vapor phase 

| = liquid phase 
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Ethyl Alcohol from 
Sulphite Waste Liquor 


A. R. KURE? 


A process used for converting the fermentable sugars in waste 
sulphite liquor into alcohol and the recovery of this aleohol is described. 
In a second process the carbon dioxide formed during the fermentation 
is used in the manufacture of magnesia insulation. 


Sprener the production of ethyl 
alcohol by fermentation of the 
sugars present in the sulphite waste 
liquor from wood pulp mills is not a 
new process, there are only three 
sulphite or cellulose alcohol plants 
operating in North America. One is 
located in the United States and two 
in Canada. However, considerable 
quantities of alcohol have been pro- 
duced during the past 50 years in 
plants located in Europe. This paper 
will describe in some detail one of 
the largest and most modern of the 
cellulose alcohol plants. 

The plant is owned and operated 
by Commercial Alcohols Limited, an 
affiliate of the Canadian International 
Paper Company, and is located ad- 
jacent to C.1.P.’s large pulp and news- 
print mill at Gatineau, Quebec. All of 
the paper mill’s recoverable su!phite 
waste liquor, totalling 5-6 million gal- 
lons per week, is utilized by the 
alcohol plant. Steam, water and elec- 
tric power are also obtained from the 
paper mill. 

The process can be divided con- 
veniently into three basic operations: 

(1) Preparation of the sulphite 

waste liquor for fermentation. 

(2) Fermentation of the sulphite 

waste liquor. 

(3) Recovery and purification of 

the alcohol produced. 


OPERATION I—Preparation of 
the Sulphite Waste Liquor for 
Fermentation 


The sulphite waste liquor as re- 


ceived from the paper mill contains 
about 10% solids, of which 15 to 20% 
are reducing sugars that are formed 
by hydrolysis of the hemicelluloses in 
the pulp ‘wood during the digester 
cooking operation. Approximately 
70% of these sugars are fermentable 
or hexose sugars, principally mannose 
and glucose. 

The liquor from the paper mill is 
received hot at about 150°F., and 
enters the first of seven wooden stor- 
age tanks connected in series. The 
liquor flows by gravity through each 
of the tanks and is finally pumped 
from the last tank into the alcohol 
plant. The seven tanks, having a total 
capacity of 900,000 gallons, or about 
one day’s production of  sulphite 
waste liquor, provide sufficient liquor 
storage to operate the distillery con- 
tinuously even though the paper mill 
is shut down for one day each week. 
The liquor remains hot, and there- 
fore sterile, during its passage through 
the storage tanks, thus minimizing the 
possibility of contamination of the 
fermentation system by wild yeasts, 
moulds or bacteria. 


Preparation of the sulphite waste 
liquor for fermentation consists essen- 
tially of cooling and neutralizing the 
liquor to the proper conditions for 
optimum fermentation efficiency. 

The liquor, which is pumped con- 
tinuously from the final storage tank, 
is cooled with water in a bank of 36 
tube-in-shell heat exchangers. The 


1 Manuscript received 


2 Assistant General Manager, Commercial Alcohols Limited, Gatineau, Que. 
This article is based on a paper which was presented at a meeting of the Pollution Control Board 
of the Province of Ontario at the Third Industrial Waste Conference, Delawana Inn, June 11, 1956. 
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final liquor temperature is controlled 
automatically to about 88°F., the op- 
timum level for fermentation. The 
resulting hot water leaving the heat 
exchangers is used as boiler feed water 
by the paper mill. The heat exchangers 
transfer about 20 to 25 million btu./hr. 
and can handle up to 40,000 gph. of 
sulphite waste liquor. 

The sulphite waste liquor flows 
from the coolers through a 3,000- 
gallon reaction vessel where the pH 
is raised from 2.5 to about 4.7 by the 
continuous addition of a 15% lime 
slurry. The reactor is divided by 
means of bafHes into three compart- 
ments, each provided with a turbine- 
type agitator. The liquor and the 
lime slurry are fed in at the top of 
the vessel, are intimately mixed, and 
flow out through the bottom. The 
flow of the lime slurry is regulated 
automatically by a continuously con- 
trolling and recording pH _ meter. 
With this system, the pH in the 
fermenters is controlled to within 
plus or minus 0.05 of a unit. 


The lime slurry is prepared and 
stored in two 6,000-gallon tanks, one 
tank being used while the other is 
recharged. Pulverized hydrated lime 
(325-mesh) is stored in a 120-ton 
capacity hopper, and is conveyed to 
the tanks through a rotary bin valve 
and screw conveyor. Recharging of 
a tank is accomplished automatically 
from a central control panel. The 
tanks are instrumented completely so 
that once the charging system is actu- 
ated the correct amounts of water and 
lime are added in the proper sequence 
and the system then shuts off. 


The hydrated lime is received in 
railway box cars and is transferred 
to the storage hopper by a pneumatic 
conveying system. A freight car hold- 
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ing 30-35 tons of lime can be emptied 
by this system in about six hours by 
one man. Approximately 20 pounds 
of hydrated lime are required for 
every 1000 gallons of sulphite waste 
liquor processed. 


OPERATION Il—Fermentation of 
the Sulphite Waste Liquor 


The cooled and neutralized liquor 
leaving the reactor flows through 40- 
mesh screens, which remove any 
coarse solids such as pulp or particles 
of unreacted lime, and then into the 
first of seven fermenters. 

Fermentation is carried on continu- 
ously under anaerobic conditions in 
seven 100,000 gallon closed wooden 
fermenter tanks connected in series. 
The fermenting liquor flows by 
gravity through the fermenters, each 
of which is equipped with a 25- 
horsepower propeller-type agitator to 
prov ide a rapid and efficient fermen- 
tation. 

A 10% slurry of yeast is fed con- 
tinuously to the first fermenter. The 
yeast, a strain of Saccharomyces cere- 
visiae that has been acclimated to a 
sulphite waste liquor environment, is 
obtained by centrifuging part of the 
completely fermented liquor or beer 
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Production of ethyl alcohol from sulphite waste liquor. 


leaving the final fermenter. The bal- 
ance of the fermented liquor, along 
with the partially clarified beer from 
the yeast separators, flows to a 100,000- 
gallon wooden beer well. The hold-up 
time in the fermenters, or the effective 
fermentation time, varies from 16 to 
21 hours, depending on the processing 
rate. Fermentation is essentially 100% 
complete at the end of the cycle, with 
about 80-85% of the fermentation 
taking place in the first two fermen- 
ters. A small amount of a 20% aque- 
ous ammonia solution is fed continu- 
ously to the first fermenter to supply 
sufficient nitrogen to maintain an ade- 
quate and healthy yeast population. 

By-product carbon dioxide gas is 
recovered from the first two fermen- 
ters and is used as a prime raw ma- 
terial in an adjacent plant where 85% 
magnesia insulation is manufactured. 

The beer, which contains about 
0.8% alcohol by volume, is pumped 
from the beer well to the distillation 
unit. 


OPERATION Ill—Recovery and 
Purification of the Alcohol 
Produced 


The alcohol is stripped from the 
beer, purified and rectified in a five 
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column distillation unit capable of 
processing up to 40,000 gph. of beer 
and producing 450 gph. of alcohol, of 
which over 85% is high purity ethyl 
alcohol (cologne spirits). Because of 
the relatively low concentration of al- 
cohol in the beer (molasses and grain 
beers average 6-9% ), steam consump- 
tion and heat utilization are important 
economic factors in the distillation 
operation. Consequently, the unit has 
been designed to obtain as high a 
thermal efficiency as possible. In addi- 
tion, features have been incorporated 
in the distillation unit to make pos- 
sible the production of cologne 
spirits. 

The five columns in the unit are 
(1) an atmospheric beer still or strip- 
ping column, (2) a purifying column, 
(3) a vacuum rectifying column, (4) 
a vacuum “heads” column, and (5) a 


vacuum “tails” column. 


The beer still is 126-in. in diameter, 
55 ft. high and operates at atmospheric 
pressure. Beer from the beer well is 
fed continuously to the beer still, after 
being pre-heated in heat exchangers 
with the hot stillage discharged from 
the base of the column. Direct low- 
pressure steam enters the bottom of 
the column and the stripped alcoholic 
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Flow splitter and yeast separators. 


vapors and other volatile materials are 
drawn from the top. Stillage, con- 
taining the non-volatile portion of the 
beer (lignosulfonates, non-fermentable 
sugars, and some dead yeast) is 
pumped continuously from the =. 
through the beer pre-heaters, and i 
discharged. 


The heat contained in the vapors 
leaving the top of the column is 
sufficient to operate the rectifying 
column and the heads and tails col- 
umns, with the only additional steam 
required being that used in the ejec- 
tors which maintain vacuums on the 
three columns. A part of the vapors 
is condensed in the calandrias of the 
heads and tails columns; the remainder 
is condensed in a vacuum steam gener- 
ator which provides sufficient steam 
(at about 15 in. Hg vacuum) to ope- 
rate the rectifying column. The con- 
densate, or high wine, from the two 
calandrias and the vacuum steam 
generator is collected in a 2,000- 
gallon drum from which it is fed to 
the purifying column. The high wine 
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alcohol. 

The purifying still is a 48-in. 
diameter hydroselection column, 60 
ft. high. The feed, preheated by the 
base discharge in tube-in-shell heat 
exchangers, enters at mid-column and 
direct steam is fed to the base. The 
feed rate is controlled automatically 
by the level in the beer still conden- 
sate drum. The heads, consisting 
mainly of volatile aldehydes, higher 
alcohols, “wood oils”, volatile sulphur 
compounds, sulphur dioxide and some 
alcohol, are condensed in a _ steam 
generator which produces a_ small 
portion of the steam that is used in 
the beer still. The non-condensibles, 
chiefly sulphur dioxide, are vented to 
the atmosphere. Part of the condensate 
is returned to the top plate of the 
column as reflux; the remainder is 
removed as a heads fraction. The 
base Sorta, containing approxi- 
mately 12% alcohol, provides the feed 
for the rectifying column. 


contains about 12 to 14% 


The purifying column functions 
both as an_ extractive distillation 


column, with water as the extracting 
agent, and as a reaction vessel. Hot 
dilution water is fed to the top plate 
of the column along with the reflux to 
effect the extractive distillation of the 
alcohol-soluble impurities. The high 
temperature and pressure in the 
column cause certain chemical reac- 
tions to take place which convert im- 
purities, normally diffiuclt to separate 
by ordinary distillation, to non-volatile 
polymerized substances that remain 
in the base discharge. These sub- 
stances are eliminated in the subse- 
quent rectification. The base draw 
from the purifying column, after 
passing through the purifying column 
feed preheaters, is cooled with water 
to the feed plate temperature of the 
vacuum rectifying column in a tube- 
in-shell heat exchanger. The cooled 
high wine is fed into the 102-in. 
diameter still at mid-column at a rate 
that is controlled by the base level in 
the purifying column. The column 
operates at a pressure of 260 mm. Hg. 
and is heated by direct steam gener- 


ated in the vacuum steam generator | 
used to condense the beer still vapors. | 


All the non-volatile material in the 
high wine feed is contained in the 
stillage that is pumped from the base 
of the column. 


The product is withdrawn from the 
upper part of the column at a strength 
of approximately 97% algohol by 
volume. The rate of draw-off is con- 
trolled automatically by the preset 
mid-column temperature. 


The overhead vapors, which con- 
tain a high concentration of volatile 
impurities, are condensed in %a tube- 
in-shell condenser. Part of this con- 
densed heads draw is returned to the 
top plate as reflux; the remainder is 
returned to the beer still condensate 
drum. 


The higher alcohols or fusel oil, 
along with some “wood oils”, ac- 
cumulate on the mid-column plates 
and are drawn off as a side cut. This 
fusel oil draw is pumped continu- 
ously to a decanter where the alcohol- 
soluble but water-immiscible oils are 
separated from the ethanol in the 
draw by dilution with water. The 
fusel oil, which is lighter than water, 
and the wood oils, which are heavier 
than water, separate as top and bot- 
tom layers respectively and are drawn 
off. The dilute alcohol middle layer 
is returned to the beer still conden- 
sate drum. 

The product from the rectifying 
column, while being a high strength 
alcohol suitable for most denatured 
grades, still contains some impurities, 
the chief one being methanol. This 
methanol originates in the digesters 
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at the paper mill, having been formed 
during the cooking process from the 
methoxyl groups in both the lignin 
and hemicelluloses in the pulp wood, 
and is present, therefore, in the sul- 
phite waste liquor produced. The 
methanol, along with small amounts 
of volatile impurities, is removed from 
the ethanol in the heads column. The 
§8-in. diameter heads column and the 
54-in. diameter tails column are each 
50 ft. high and operate at about 200 
mm. Hg. vacuum. They both are 
heated indirectly by condensing some 
of the beer still vapors in a calandria 
in the base of each column. 

The product from the rectifying 
column is fed into the upper part of 
the heads column. The methanol i 
the feed, having a lower boiling point 
than ethanol, is concentrated in the 
overhead vapors along with the re- 
maining volatile impurities such as 
aldehydes, and is removed in the 
heads draw as an anhydrous methanol- 
ethanol mixture. This draw is con- 
trolled automatically by the preset 
temperature at the top of the column. 
Part of the heads draw is returned to 
the top plate as reflux while the re- 
mainder is stored for use as a denatur- 
ant or for sale as a solvent. 

The base discharge from the heads 
column is at about the same alcohol 
strength as the feed to the column 
but now contains no methanol. The 
heads column product is fed to the 
tails column for the final purification 
step, the rate of feed being controlled 
automatically by the base level in the 
heads column. 

The feed enters the final still below 
mid-column and the cologne spirits 
product is withdrawn under auto- 
matic control. Any last traces of 
volatile impurities that might be pre- 
sent in the feed to the column are 
concentrated in the overhead vapors 
which are condensed in a tube-in-shell 
condenser. Part of the condensate is 
returned to the top plate as reflux 
and the remainder is sent to the beer 
still condensate drum for recycling. 
Any remaining higher boiling impuri- 
ties, such as higher alcohols and cer- 
tain esters, are concentrated in the 
base and are removed as a tails draw 
and stored for use as a denatured 
grade product. The draw rate is con- 
trolled automatically by the base level. 

The methanol draw, the alcohol 
draw for denatured grade product, 
and the cologne spirits draw are col- 
lected in separate closed receivers 
in accordance with the Government 
Excise regulations. From the receiv- 
ers, the various products, after having 
been approved for quality by the 
control laboratory, are w eighed and 
tested in the presence of Government 


Excise Officers and then pumped to 
storage. 


Process and Quality Control 

In order to maintain continuous 
production of a high quality ethyl 
alcohol, the distillery is instrumented 
as completely and ‘thoroughly as is 
practical. Nearly all the instruments 
are operated pneumatically and are 
located at three control stations: one 
station in the sulphite waste liquor 
cooling and neutralizing area, another 
in the fermentation area and a third 
in the distillation area. By means of 
these instruments the entire distillery 
is started up, operated and shut down 
by only three operators per shift. 

While the continuous operation 
and extensive instrumentation insures 
relatively uniform product quality, 
the intermittent batch operation of 
the digesters at the mill tends toward 
variations in the sulphite waste liquor 
composition from cook to cook. In 
addition, changes are made in the 
grade of pulp produced, the type of 


wood used, and the cooking condi- 
tions. These variables affect the per- 
cent of solids, sugars, free and com- 
bined sulphur dioxide, methanol, and 
calcium sulphate present in the liquor. 

To minimize the effect of these 
variations, the quality control labo- 
ratory carries out routine chemical 
and phy sical analyses on the various 
process streams in the alcohol plant. 
The results of these analyses are used 
to determine the optimum operating 
conditions for the liquor being pro- 
cessed. Yeast cell count and viability 
tests are performed regularly to in- 
sure maximum fermentation efficiency. 
Routine bacteriological tests are also 
carried out on the liquor at various 
stages in the process to guard against 
bacterial contamination of the fer- 
mentation system. 


In addition to these analyses each 
receiver of the final ethyl alcohol 
product is tested for product quality 
before being put into storage. An 
alcohol that does not meet the Com- 





Sulphite waste liquor coolers. 
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pany specifications for cologne spirits, 
which are more rigid than both the 
U.S. and British Pharmacopoeia 
specifications, is put into storage for 
denatured grade alcohol. 

Depending on the sugar content of 
the sulphite waste liquor, a total of 

2.5 to 3 million gallons of alcohol can 
be produced annually, of which about 
85% is high quality cologne spirits. 


Production of 85% Magnesia 
Insulation 

About 2.5-3 million board feet of 
85% magnesia insulation are produced 
annually for Canadian Johns-Manville 
in an adjacent plant. This plant utilizes 
some of the 150-190 million cu. ft. 
of by-product carbon dioxide pro- 
duced annually by the distillery. It is 
interesting to note that about 30% of 
the weight of the final product is pro- 
vided by this carbon dioxide. 

In this process the fermenter gases 
from the first two fermenters, where 
about 80-85% of the fermentation 
takes place, are pumped directly to 
a 3,500-gallon continuous carbonator 
in the magnesia plant. The carbon- 
ator is charged continuously with 
water, a solution of magnesium bicar- 
bonate, and finely pulverized magne- 
sium hydroxide. The bicarbonate and 
hy droxide react to form magnesium 
normal carbonate crystals, and the 
fermenter gases, containing about 
95% carbon dioxide, are dispersed in 
the carbonator by a high speed agi- 
tator to maintain the required bicar- 
bonate concentration. The dilute sus- 
pension of magnesium normal car- 
bonate crystals overflows continuous- 
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ly to an Oliver vacuum filter where 
the crystals are concentrated and put 


into storage as a 20% slurry. The 
filtrate, which contains magnesium 
bicarbonate, is returned to the car- 


bonator. To control the carbonation 
rate, a constant bicarbonate concen- 
tration is maintained in the carbon- 
ator by adjusting the feed rates of 
the hy droxide and make- up water. 


To make the molding slurry, the 
20% suspension of magnesium normal 
carbonate crystals is blended with 
asbestos fibre in 400-gallon batches. 
The asbestos is added as a temper- 
ature-resistant binder to give added 
strength to the final product. 


Just prior to molding, the carbon- 
ate-fibre slurry is diluted with hot 
water which partially converts the 
magnesium normal carbonate crystals 
to magnesium basic carbonate. This 
reaction takes from eight to ten min- 
utes, after which the slurry is molded 
into insulation. 


Two molders of the reciprocating 
filter press type, designed by Johns- 
Manville, and capable of molding 
three-foot lengths of pipe covering 
or block insulation of various sizes 
and thicknesses, are required for pre- 
sent production. A metered charge of 
dilute magnesium basic carbonate and 
asbestos fibre is dumped into the mold 
box which rests on one of two reci- 
procating pallets covered by wire 
mesh. A ram, to which is attached a 
perforated steel molding head also 
covered with wire mesh, descends 
into the mold box and the excess 
water is pressed out of the slurry. 


The plastic-like residue of magnesium 
basic carbonate and asbestos fibre 
takes the shape of the space between 
the pallet, the molding head and the 
mold box sides. 


After a piece has been molded, the 
ram and the mold box rise, the pallet 
with the freshly molded piece swings 
from under the mold box and the 
empty pallet moves into position. 
Then the piece is removed from the 
pallet, placed on a curing table, and 
later placed in a drying car. While 
this is taking place the mold _ box 
descends, a fresh charge is metered 
into the box, and the molding cycle 
is repeated. 


The molders are hydraulically ope- 
rated, and the molding sequence is 
controlled automatically by a series 
of electrically driven adjustable cams 
that open and close the various valves 
in the hydraulic system, either direct- 
ly or by triggering solenoid valves. 
Each molder is operated independent- 
ly through its individual control sys- 
tem. 


The molded insulation is dried for 
about 24 hours in forced-air tunnel- 
driers at about 300-340°F. The dried 
blocks are removed from the drier 
cars, inspected and packed in car- 
tons for shipment. The dried pipe 
covering, made in half sections, must 
be trimmed so that the two pieces 
provide a uniform fit on standard steel 
pipe. Then the pieces are inspected, 
made into whole sections by w rapping 
with cotton cloth, and packed i in Car- 
tons for shipment. 


¥ -* -»* 
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Underground Storage Created 


in Salt Beds at Sun’s 


Sarnia Refinery 


G. R. HENDERSON? and P. F. DOUGHERTY* 


One method of storing liquified petroleum gases in underground 


salt bed caverns is described. 
the cavern are discussed. 


HE chemical industries at Sarnia 
"Laan to a considerable extent 
on light hydrocarbons obtained from 
the local refineries as feedstock for 
their operations. The most desirable 
hydrocarbons are the C, (ethylene- 
ethane), Cs (propylene- propane), and 
C, (buty rated butane) fractions and 
the synthetic rubber plant has re- 
quired an ever increasing supply 
for the Cy (butane, butylene mixture ) 
since it commenced operation. Imme- 
diately after the termination of World 
War Il the refining capacity in the 
area was not sufficient to meet the C, 
requirements of the customer and for 
many years this material was imported 
by railroad tank car from the Sun Oil 
Co. refinery at Toledo. 

The transportation of large volumes 
of butane and propane by tank car, 
is expensive and, in addition, the 
shipping and receiving companies have 
to invest large capital expenditures 
in loading racks and pressure tankage. 

One of the obvious solutions to the 
transportation problem was the ship- 
ment of the C, fraction from Toledo 
to Sarnia by pipe line. When the 
Interstate Commerce Commission de- 
clared the majority of pressure tank 
cars obsolete in 1950, Sun Oil elected 
to install a pipe line between Toledo 
and Sarnia for the shipment of butanes. 
The Company intended, eventually, 
to build a refinery and then the pipe 


Some of the factors involved in preparing 


line would also be available for crude 
deliveries. The initial operation 
through the pipe line was a simple one 
product movement, and the line itself 
supplied some 40,000 bbl. of storage 
capacity. However, with the opening 
of the Sarnia Refinery it was neces- 
sary to make further shipments by 
“slugging” some 15,000 to 20,000 bbl. 
of C, at a time between crude deliver- 
ies. By this method of transportation 
the C, arrives at Sarnia as a dirty 
“BB” fraction because of the con- 
tamination with crude oil at either end. 

Conventional spheres or spheroids 
for storing the “BB” fraction pre- 
sented numerous problems due to 
initial cost and method of operation. 
However, the experience that had 
been gained with underground salt 
cavity storage of LPG in Texas about 
1950 indicated that this method of 
providing adequate storage at Sarnia 
was a real possibility. Local geology, 
obained from drilling operations in 
the district indicated the presence of 
continuous beds of salt underlying 
the refinery area and the original core 
drilling on the first well indicated 





Type of Storage 


High- pressure above ground storage (for Propane) 


Mined Storage (for Propane- -Butane) 


Low-Pressure-above ground - Sphereoid- -refrigerated (for Butane) 


Salt Cavern (for Propane-Butane) . 


1 Manuscript received April 18, 1957. 


2 Manager, Sarnia Refinery, Sun Oil Co. Limited, Sarnia, Ont. 
3 Process Engineer, Manufacturing Department, Sun Oil Co. Limited, Sarnia, Ont. 


This article is based on a paper which was presented at the March 11, 


1957 meeting of the 


Chemical Engineering Division of The Chemical Institute of Canada. 
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the existence of two salt beds which 
could be used. The first, encountered 
at about 2,000 ft., was 200 ft. thick 
and the second, which started about 
2500 ft. had a thickness of 137 ft. and 
was comparatively free from in- 
solubles. 

The original geological data ob- 
tained indicated that storage cavities 
washed out of the salt would be free 
from leakage. Subsequent operating 
data has confirmed this original as- 
sumption. 

The underground storage at Sarnia 
was installed for the specific purpose 
of transferring C, from a pipe line 
shipment to the synthetic rubber 
plant. This contaminated material 
which is stored in large quantities as 
dirty “BB” is cleaned by fractionation 
at the refinery and delivered to the 
customer at a rate of about 1,000 bbl. 
a day. The cavities are also used to 
receive the returned saturated butane 
portions of Sarnia refinery’s produc- 
tion. Large storage capacity provides 
another advantage to the refinery 
because LPG can be retained when 
off season demands create non-pro- 
fitable prices for the product. A com- 
parison of the costs of storing light 
hydrocarbons by various methods is 
as follows: 


Cost per Bbl. Dollars 


20.00 - 30.00 
7.00 - 15.00 
12.00 - 20.00 
0.75 - 3.50 


In addition to the above variation 
in cost it has been found that the 
operation and maintenance of the salt 
cavities is much less than on the con- 
ventional type of equipment. To in- 
dicate the trend toward underground 
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storage the petroleum industry has 
developed about 15,000,000 bbl. of 
storage in salt strata during the past 
five vears. 

In determining the feasibility of 
salt cavern storage the following fac- 
tors were considered: 

(1) the depth of the salt beds. 

(2) The availability and cost of 
fresh water for washing, and 
means of disposing of the brine 
resulting from washing. 

(3) the possibility of materials asso- 
ciated with the salt having a 
harmful effect on the material 
being stored. 

A preliminary study indicated that 
all the above factors presented no 
problem. The salt beds were at a 
reasonable depth, no problems existed 
regarding the supply of fresh water 
or brine disposal, and the salt beds 
contained no materials which would 
have a deleterious effect on the ma- 
terials which were being stored. 

As a result, five caverns have been 
drilled and washed. At the present 
time the total capacity of these wells 
is 450,000 bbl. Present plans include 
the washing of two of these wells to 
approximately 200,000 bbl. each, rais- 
ing total capacity to something in the 
order of 650,000 bbl. 

When it was decided to create 
storage cavities at Sarnia the follow- 
ing had to be determined: 

1—Spacing pattern 

2—Drilling and Casing program 

3—Washing program 

4—Methods for testing the tubing 

5—Methods for strapping the 

caverns 

6—Brine pit operation 


Spacing Pattern 
The spacing pattern of the wells at 

Sarnia was based on reducing the 

chance of intercommunication — be- 

tween the finished caverns as a result 

of washing or accidental overpressur- 

ing, and also on the economy of the 

interconnecting piping between the 

caverns. 

The following assumptions govern- 

ed cavern spacing. 

First —The ultimate capacity of each 
cavern would be 200,000 bbl. 

Second—The effective thickness of the 
salt bed which was washed 
would be 100 ft. (The over- 
all thickness of the lower salt 
bed is approximately 127 ft.). 

Third —The diameter, calculated from 
the two previous points, is 
approximately 120 ft. This 
was doubled on the assump- 
tion that a cavity could be 
washed to one side of the 
piping. 


Inasmuch as a minimum space of 
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350 ft. between the caverns was de- 
sired, the bore holes through which 
the caverns were to be washed were 
drilled 600 ft. apart. 


Drilling and Casing Program 

A 173 in. hole was drilled to a 
depth of approximately 200 ft. into 
which a 12} in. casing was set and 
concreted into the surface of the 
ground (see Figure 1). Using this 
casing an 11 in. hole was drilled with 
fresh water mud to a depth of 1,425 
ft. At this point salt-saturated mud 
was used and a 6 in. hole was con- 
tinued completely through the salt 
formation. This hole was reamed to 
11 in. and an 88 in. casing was set 
from the top of the desired cavity 
and cemented in place. 

After drilling out the casing plug, 
hydrostatic pressure was placed on 
the casing seat and the open hole be- 
low, to a maximum pressure of 0.8 psi. 
per ft. of depth of the hole. This is 
the pressure at the bottom of the hole. 

The enlargement of the 6 in. hole 
to 11 in. was then continued to the 
desired total depth. A 7 in. floating 
string was run to the depth required 
to afford, by means of a kerosene seal, 
protection to the roof of the cavern 
during the washing operation. Then 
4} in. Hydril tubing was run to with- 
in approximately 3 ft. of the bottom 
of the drill hole. 


Washing Program 


The factors which normally must 
be considered in determining the 


optimum rate at which a cavern 
should be washed are: 
1—Availability and cost of fresh 


water. 

2—Problem of brine disposal. 

3—Thickness of the salt bed being 

washed. 

4—Nature of the insoluble material 

in the salt bed. 

5—Time limitations. 

At Sarnia there was no problems 
created by the suppy or cost of the 
fresh water, brine disposal, or the 
nature of the insoluble materials in 
the salt beds. For this reason, the 
prime consideration was the time ele- 
ment, inasmuch as the caverns were 
required for refinery operations at a 
specified date. 

The first 5,000 bbl. of cavity were 
made by injecting fresh water through 
the 44 in. Hydril tubing, and the re- 
sulting brine removed through the 
annulus between the 43 in. tubing 
and the 7 in. floating string. Although 
this was the reverse of normal pro- 
cedure and an inefficient method of 
washing, it theoretically provides 
space for the insolubles which may 
settle out as washing proceeds. The 


main washing was performed by ad- 
mitting water through the annulus 
and withdrawing the brine through 
the tubing. 

During the preliminary washing, 
sufficient kerosene was injected to fill 
the annulus between the floating 7 in. 
string and the 8 in. inner casing to 
give some measure of protection to 
the roof. After the preliminary wash- 
ing the cavern was pressure tested to 
600 psig. at the surface and then the 
flow was reversed to the normal pro- 
cedure. 


Reasons for Testing the Tubing 
Past experience indicated the need 
for establishing the tightness of the 
tubing before the cavern is placed in 
use. 
|—Leakage of the tubing can result 
in liquefied gas entering the tub- 
ing from the annulus and large 
volumes might be released to 
the atmosphere. 
—There is a further possibility that 
a leak could indicate complete 
breakage of the tubing during 
a period of liquefied gas injection. 
As a matter of safety, such leak- 
age would require the cavern to 
be emptied and the tubing re- 
paired. 


Methods of Strapping Caverns 

In view of the method employed 
in washing the cavern and the relative 
absence of large stringers, it can be 
assumed that the shape of the cavern 
approximates a cylinder. Thus, a given 
volume-height relationship will fix the 
approximate diameter. This volume- 
height relationship can be obtained by 
measuring the difference in shut-in 
annulus pressures of the liquefied light 
hydrocarbons at the surface after 
successive additions of known volumes 
to the cavern. The rise in pressure 
together with the calculated base 
pressure required on the annulus to 
balance the brine column from the 
bottom of the casing seat permits cal- 
culation of an equivalent diameter at 
the various test points. 

The maximum diameter of the 
cavern calculated in this manner is 
useful since it provides a check on 
the validity of the assumed well- 
spacing pattern. 


Brine Pit Operation 

To reduce the hazard resulting from 
a release of hydrocarbons through the 
brine line to the brine pit, a flare is 
kept in operation at all times at the 
outlet of the brine line. A sudden 
release of hydrocarbons may result 
from a break in the tubing or from 
accidental overfilling of the cavern. 
The flare ignites any hydrocarbons 
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thus vented for obvious safety reasons 
and the value of this procedure has 
bec:: confirmed by operating experi- 
enc: 

When the cavity overflows, hydro- 
carbons reach the bottom of the 44 
in. tubing and the seal is lost on the 
column of brine which drains out of 
the cubing. The result is an equaliza- 
tion of pressure on either side of the 
tubing and an increase in volume of 
the hydrocarbon in the cavity be- 
cause of the reduction in pressure. 
The volume increase is very appreci- 
able and has to be taken out of the 
well before normal conditions can be 
re-established. This increase in vol- 
ume has to be burned at the flare 
stack and will amount to several hun- 
dred barrels. To say the least this is 
an expensive procedure. 

From the foregoing it will be real- 
ized that it is important to maintain 
accurate accounting records so that 
these occurrences can be kept to an 
absolute minimum or entirely  eli- 
minated, From experience gained dur- 
ing the operation of these cavities it 
is well-known that accurate instru- 
mentation of receipts and deliveries is 
a very essential factor in maintaining 
satisfactory performance. 

The static pressure caused by the 
difference in specific gravity between 
the brine and the C, is approximately 
600 Ib. at the Sarnia cavities. There- 
fore, it is necessary to have a pump 
with about 1,000 ib. discharge pres- 
sure to overcome the static head and 
line loss to permit filling rates equiva- 
lent to pipe line deliveries. In receiv- 
ing pipe line deliveries, the cut be- 
tween crude oil and C, is made by a 
mid-gravity point, pumping is stopped 
until valve changes can be made at 
the crude tanks and the well head. A 
similar procedure is adopted at the 
tail-end of the ee Some crude 
is mixed with the C, and the material 
is cleaned by simple distillation as it is 
removed from the dirty BB well and 
before it is delivered to the customer. 

At times customer requirement may 
exceed the c capacity of the distillation 
column and it was found that a clean 
“heart cut” of about 60% of the ship- 
ment could be removed to a clean 
well. This procedure permits increas- 
ing delivery rates to the customer 
and eliminates the cost of re-running 
operations. 

The first time that the wells are 
filled and emptied, some Cy, is lost in 
permanent storage. It is believed that 
this occurs because material may be 
washed above the desired cavity roof 
and, therefore, above the bottom of 
the annulus pipe. Experience over a 
period of several years indicates yp 
such permanent storage loss is of < 
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Figure 1—Cross-section view of one of the salt caverns at Sun Oil’s Sarnia 


refinery show the displacement of casing and tubing. 


Washed to approximately 


2,600 ft., all five wells at this location provide a total storage capacity for about 


450,000 bbl. of light hydrocarbons. 


relatively small volume. In one of the 
more recent wells, which was pro- 
vided with an oil cap to prevent the 
water reaching the top of the cavity, 
the permanent inventory after the 
first filling and emptying was found 
to be so small a TT that it could 
be considered a gauging error. 

Salt cavities have been used success- 
fully for both C, and gasoline. Ex- 
haustive tests and actual operations 
by Sun Oil indicate that there is no 


detrimental effect on the stored 
material. 
Underground storage at the Sarnia 


refinery has satisfied all the require- 
ments for which it designed. 
However, successful operation is only 
achieved by experienced operators 


was 


and good maintenance, especially on 
the vital instrumentation, which pro- 
vides the only record of the material 
transfers. : 

The movement of hydrocarbons in 
and out of the cavities and involving 
pipe line deliveries was without pre- 
cedent and local procedure had to be 
developed by trial and error. 

The constant movement of light 
hydrocarbons, without the advantage 
of conventional gauging procedure, 
has been a trying experience, but the 
advantages of lower initial cost and 
the safety of underground storage, 
have far outw eighed the periodic 
complicated problem of determining 
current inventory. 
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